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1  INTRODUCTION 


The  United  States  Department  of  Energy  (DOE)  in  1985  established  the 
Subsurface  Science  Program.  The  primary  purpose  of  this  program  is  to  assist 
doe’s  efforts  to  clean  up  wastes  that  have  been  disposed  of  underground 
and  that  have  moved  from  surface  sites  into  the  underground.  While  the 
mining  and  petroleum  companies  have  carried  out  extensive  operations  in 
the  subsurface,  the  physics,  chemistry,  and  particularly,  the  biology  of  the 
subsurface  environment  are  largely  unknown. 

In  Section  2,  we  present  our  overall  view  of  the  program,  its  major  goals, 
and  outstanding  issues.  In  general,  this  is  a  high-quality  program,  and  de¬ 
serves  continuing  and  increasing  support.  The  problem  of  understanding  the 
subsurface  environment  in  such  a  way  as  to  make  bioremediation  a  useful  al¬ 
ternative  is  highly  challenging.  The  subsurface  environment  is  characterized 
by  great  natural  physical  and  biological  heterogeneity.  This  heterogeneity  is 
discussed  in  Sections  2,  3,  and  6.  The  complexity  of  natural  heterogeneity  is  a 
theme  running  throughout  the  report.  Methods  to  describe  the  heterogeneity 
in  a  useful  way  must  be  developed. 

The  questions  raised  in  Section  2  attempt  to  frame  the  goals  of  a  sub¬ 
surface  program.  These  goals  should  include  understanding: 

•  subsurface  movement 

•  subsurface  microbial  structure 
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•  in  situ  behavior  of  subsurface  microbes 


•  origins  of  subsurface  microbial  organisms. 

In  pursuing  these  objectives,  there  should  be  a  balance  between  tech¬ 
nique  development  and  field  studies. 

Of  the  various  issues  facing  the  Subsurface  Program,  we  have  selected 
one — the  transport  of  bacteria  underground — for  a  more  detailed  technical 
treatment.  Section  3  deals  with  the  physics  of  attachment  of  bacteria  to 
walls,  while  Section  4  considers  the  problems  in  understanding  transport  of 
bacteria  through  porous  media.  In  selecting  these  problems,  we  recognize 
that  we  do  not  yet  know  how  bacteria  function  in  the  subsurface,  whether  in 
situ  or  planktonic  (free). 

Section  5  provides  a  brief  introduction  to  bioprocessing  of  metals.  The 
subject  has  a  vast  literature,  and  therefore  we  only  present  an  overview. 

In  Section  6,  we  consider  the  issue  of  transmitting  genetic  information 
to  resident  bacteria.  We  also  discuss  the  possibility  of  using  viruses  (phage) 
and/or  liposomes  to  provide  the  genetic  information  to  resident  bacteria  so 
as  to  give  them  the  desired  properties.  Section  2  comments  on  the  validity 
of  sampling  the  subsurface. 

Section  8  departs  from  the  main  theme  of  the  report.  In  this  section, 
we  briefly  consider  the  difficult  issue  of  nuclear  waste  disposal.  The  reason 
for  doing  so  is  to  emphasize  that  alternatives  to  borosilicate  glass  exist,  and 
that  these  alternatives  should  not  be  forgotten,  particularly  because  long¬ 
term  disposal  has  been  greatly  delayed. 
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In  Section  9,  we  comment  on  the  new  Environmental  Molecular  Sci¬ 
ence  Laboratory  (EMSL).  We  believe  that  this  new  laboratory  can  be  one 
of  doe’s  centers  supporting  subsurface  science.  It  is  clearly  important  that 
the  EMSL  have  microbiological/molecular  biological/biochemical  expertise 
in  residence.  However,  we  are  concerned  that  the  Laboratory’s  mission  is 
not  well  defined. 

In  the  final  section,  we  briefiy  consider  the  long-term  future  of  bioreme¬ 
diation.  A  science  in  its  infancy  probably  should  not  go  into  the  field,  but  the 
future  of  bioremediation  as  an  option  depends  on  field  investigation.  There 
thus  needs  to  be  the  balance  discussed  above  between  field  work,  laboratory 
investigation,  and  methods  development. 

The  report  does  not  pretend  to  be  a  comprehensive  overview  of  the  entire 
Subsurface  Program.  For  example,  we  say  nothing  about  ongoing  work  on 
issues  related  to  mixed  contaminants.  What  we  have  tried  to  do  is  provide 
an  introduction  to  some  important  elements  of  the  program,  and  to  discuss 
some  ways  in  which  the  program  could  be  improved. 
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2  SUBSURFACE  SCIENCE  PROGRAM 
OVERVIEW 


In  general,  the  Department  of  Energy  (DOE)  Subsurface  Science  Pro¬ 
gram  (SSP)  is  an  excellent  one  of  high  scientific  quality.  It  is  unique  in  the 
degree  to  which  the  various  disciplines  have  been  brought  together,  and  in 
the  integration  of  field  and  laboratory  programs.  In  many  ways,  it  is  an  out¬ 
standing  paradigm  for  future  programs  that  must  integrate  basic  and  applied 
sciences  from  several  disciplines  in  the  service  of  highly  practical  goals. 

The  critical  issue  for  the  future  of  the  program  is  to  have  a  clearly 
articulated  goal.  This  goal  can  then  be  used  to  judge  the  relevance  of  var¬ 
ious  scientific  activities.  Since  the  subsurface,  and  in  particular  subsurface 
ecology,  is  such  an  unexplored  frontier,  there  may  be  a  tendency  to  pursue 
approaches  of  high  scientific  interest  that  do  not  necessarily  contribute  to  the 
overall  goal  of  the  program.  Our  understanding  of  the  long-term  practical 
goal  of  the  program  is  to  enable  and  enhance  the  capability  of  the  United 
States  to  diagnose,  understand,  and  carry  out  effective  and  economically  fea¬ 
sible  subsurface  remediation.  In  furthering  this  goal,  there  is  a  real  need 
to  answer  a  number  of  open  scientific  questions.  While  we  do  not  claim 
completeness,  the  following  questions  appear  to  have  major  significance  for 
reaching  the  overall  practical  goal. 

1.  How  can  we  acquire  a  fundamental  (predictive)  understanding  of  the 
movement  of  fluids,  chemicals,  and  colloids  (including  biological)  in 
complex  subsurface  media?  It  is  clear  that  our  present  ability  to  pre¬ 
dict  subsurface  movement  is  minimal  at  best.  A  major  part  of  this  is 
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physics  and  chemistry,  but  a  significant  section  is  biological.  The  bio¬ 
logical  component  of  subsurface  flow  is  largely  unexplored.  Elsewhere 
in  this  report  we  present  a  technical  discussion  of  the  movement  of 
microorganisms  in  the  subsurface  environment. 

2.  What  microorganisms  are  present  in  the  subsurface?  An  understanding 
of  the  microbiology  census  of  the  subsurface  would  seem  to  be  a  prere¬ 
quisite  to  the  understanding  of  what  microorganisms  do,  and  are  capa¬ 
ble  of  doing,  to  the  subsurface  environment.  Since  the  rock  structures  of 
the  subsurface  environment  are  highly  heterogeneous  and  anisotropic, 
there  are  difficult  problems  of  how  to  construct  a  statistically  sensible 
sampling  strategy.  Conventional  statistics  are  probably  inadequate  to 
characterize  the  subsurface  environment  in  a  meaningful  way.  Explicit 
account  must  be  taken  both  of  the  anisotropy  of  the  subsurface  and  of 
the  fact  that  variations  at  virtually  all  scales  are  present. 

3.  How  do  known  microorganisms,  much  less  unknown  species,  behave  in 
the  subsurface  environment?  What  chemical  transformations  do  they 
accomplish,  under  what  conditions,  and  what  past  modifications  of  the 
environment  are  attributable  to  these  bacterial  populations?  All  of 
these  questions  concern  scientific  uncertainties  that  are  central  to  the 
program. 

4.  What  are  the  origins  of  the  subsurface  microorganisms?  Have  they 
developed  in  place,  or  have  they  been  carried  from  the  near-surface 
layer  to  the  deep  subsurface?  To  what  extent  have  they  evolved  in 
their  present  state?  What  do  these  subsurface  microorganisms  tell  us 
about  the  origin  of  life? 
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Given  that  a  large  part  of  the  significant  scientific  uncertainty  concerns 
the  biological  component  of  the  subsurface,  and  given  that  the  success  of 
the  program  depends  on  reducing  this  uncertainty,  it  makes  sense  that  a 
significant  effort  be  devoted  to  the  biological  components  of  the  program. 
In  particular,  a  much  more  concerted  effort  should  be  developed  to  adopt, 
develop,  and  utilize  modern  molecular  biological  methods  for  the  subsurface 
program.  For  example,  methods  such  as  those  developed  by  David  White 
of  the  University  of  Tennessee  to  characterize  the  signature  lipid  compounds 
made  by  various  organisms,  and  complementary  DNA-based  methods,  could 
significantly  increase  the  sensitivity  of  detection  and  the  density  of  informa¬ 
tion  obtained  about  subsurface  organisms.  The  polymerase  chain  reaction 
method  (PCR),  for  instance,  has  been  shown  to  be  capable  of  amplifying, 
detecting,  and  identifying  a  single  DNA  molecule,  which  translates  to  a  single 
bacterial  cell  if  its  DNA  is  accessible.  PCR  is  routinely  used  on  marine  mi¬ 
crobial  samples.  It  has  also  been  used  on  in  situ  biological  samples  of  various 
kinds,  such  as  fixed  tissue  sections.  It  should  be  possible,  with  little  overall 
effort,  to  develop  the  techniques  for  using  PCR  effectively  and  inexpensively 
on  subsurface  samples.  This  would  change  dramatically  what  could  be  ob¬ 
tained  from  the  field  samples.  While  the  program  is  already  devoting  some 
effort  to  the  development  of  these  techniques  for  using  DNA-based  methods 
(DNA  probe  methods,  for  example),  they  should  be  a  priority  of  the  program. 
The  methods  are  increasingly  versatile  and  powerful  for  characterizing  bio¬ 
logical  organisms.  Only  with  such  methods  can  the  field  programs  gain  their 
full  effectiveness  and  cost  efficiency  in  knowledge  gained  for  dollars  spent. 

Thus,  there  should  be  a  better  balance  between  the  development  and 
use  of  new  laboratory  techniques  adaptable  to  the  field  and  the  field  stud¬ 
ies  component  of  the  program.  It  will  also  be  necessary  to  perform  more 
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biochemical  and  genetic  characterization  of  subsurface  organisms  as  part  of 
the  program  in  order  fully  to  realize  the  potential  gain  of  useful  knowledge. 
These  studies  should  include  some  selective  sequencing  of  novel  organisms’ 
genomes  and  some  biochemical  physiology  of  the  bacteria,  to  include  studies 
of  the  control  mechanisms  for  gene  expression  operative  in  the  subsurface. 
This  latter  is  too  ambitious  a  goal  for  the  program  to  approach  on  a  broad 
front,  but  selective  efforts  and  a  concerted  effort  to  open  the  lines  of  commu¬ 
nication  and  collaboration  of  the  program  with  other  microbiologists  working 
on  such  problems  should  be  a  priority.  The  geologic  component  of  the  pro¬ 
gram  and  the  origin  of  subsurface  bacteria  effort  are  exciting  scientifically, 
generally  well  conceived,  and  should  lead  to  major  modifications  of  our  in¬ 
sight  into  the  origin  of  subsurface  microorganisms.  This  insight  is  likely  to 
have  a  significant  practical  outcome. 
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3  ATTACHMENT  OF  BACTERIA  TO  WALLS 


3.1  Introduction 


Toxic  waste  that  has  had  long  times  to  percolate  through  the  subsurface 
presents  major  problems  for  any  type  of  clean-up  effort.  In  bioremediation, 
the  transport  of  bacteria  to  the  subsurface  sites  of  the  waste  and  insurance 
of  bacterial  activity  once  on-site  are  both  major  issues.  The  nature  of  bac¬ 
terial  adhesion  to  surfaces  is  a  major  factor  in  both  of  these  problems,  and 
adhesive  properties  thus  will  be  an  important  criterion  in  selecting  and/or 
bio-engineering  bacteria  for  subsurface  waste  degradation. 

The  evaluation  of  how  bacterial  adhesive  properties  will  affect  func¬ 
tion  in  biodegradation  is  a  difficult  problem  due  to  the  many  length  scales 
of  heterogeneity  involved.  Bacterial  surfaces  themselves  have  complex  and 
heterogeneous  structures  [1]  [2],  so  that  the  simple  physical  model  [1]  [3] 
generally  used  to  describe  their  temporal  adhesive  behavior  has  limited  pre¬ 
dictive  capability.  In  addition,  the  process  of  bacterial  adhesion  is  likely  to 
be  temporally  heterogeneous,  as  the  bacteria  respond  biologically  to  changes 
in  their  environment.  Finally,  the  mineral  walls  with  which  the  bacteria  in¬ 
teract  below  ground  are  highly  heterogeneous,  both  chemically  and  in  terms 
of  roughness,  on  lengths  ranging  from  the  atomic  to  the  macroscopic. 

The  DOE  Subsurface  Science  Program  is  addressing  many  of  these  is¬ 
sues  through  a  broad-based  program  including  microbiological  analyses  of 
indigenous  bacteria,  laboratory  studies  of  adhesion,  field  studies  of  geologi- 
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cal  sites,  and  field  and  laboratory  studies  of  transport  in  mixed  media.  The 
results  of  all  of  these  efforts  need  to  be  coordinated  in  order  to  establish  a 
successful  bioremediation  effort. 

In  the  following,  we  will  review  some  of  the  basic  factors  determining 
cell  wall  adhesion,  and  discuss  what  must  be  known  about  these  factors 
for  prediction  of  transport  processes  in  porous  media.  We  end  with  some 
suggestions  about  how  the  existing  effort  to  develop  successful  bioremediation 
strategies  can  be  strengthened. 

3.2  Physical  Basis  of  Cell- Wall  Interactions 


Attachment  forces  between  a  bacterium  and  a  wall  are  difficult  to  calcu¬ 
late,  and  perhaps  even  estimate  with  confidence,  from  first  principles.  They 
have  several  origins:  electrostatic  interactions,  van  der  Waals  forces,  polymer 
bridges,  through  covalent  chemical  bonding;  and  other  forms  of  contact  chem¬ 
ical  bonding.  The  attachment  can  be  reversible  (the  detached  bacterium  and 
the  wall  do  not  differ  from  pre-detachment  ones)  or  irreversible  (the  detached 
bacterium  leaves  behind  a  polymer  “footprint”  on  the  wall).  Present  theo¬ 
retical  estimates  for  the  probability  and  strength  of  bacterial  attachment  to 
walls  must  rely  on  gross  oversimplifications,  which  can  make  extrapolations 
to  new  wall  materials  and  conditions  unreliable  until  measured. 


3.2.1  Electrostatic  Interactions 


In  water,  bacteria  generally  have  a  negative  surface  charge  and  surface 
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potential  (VJ,).  Similarly,  mineral  surfaces  in  contact  with  water  (e.g.,  glass, 
clay)  also  acquire  negatively  charged  surfaces  and  potentials  (Vj).  The  elec¬ 
tric  field  from  both  surfaces  is  screened  by  the  positive  and  negative  ions 
of  the  water  (Debye  screening),  especially  if  it  contains  additional  dissolved 
electrolytes.  For  weak  potentials  {eVb,a  hr) 

-  A-V  =  0.  (3  -  1) 


The  Debye  shielding  length 


(  hTe 


1/2 


where  is  the  number  density  of  dissolved  ion  pairs  and  e  is  the  dielectric 
constant  of  water.  For  distilled  water  (n^  ~  10^^  cm“^)A  is  about  the  same 
as  the  radii  of  typical  bacteria;  it  is  very  much  less  than  those  radii  in  the 
cases  of  interest  here  where  the  electrolyte  concentration  is  very  much  higher. 


The  “standard”  model  assumes  that  as  a  bacterium  in  water  approaches 
a  wall-water  interface  the  potentials  of  both  surfaces  remain  what  they  were 
when  the  bacterium  was  infinitely  far  away.  Then,  for  an  assumed  spherical 
bacterium  of  radius  R{^  A)  with  a  distance  of  closest  approach  h  between 
its  surface  and  flat  wail  (Figure  3-1),  the  Coulomb  interaction  energy  is  [3] 

[4]  [5] 

E.{h)  =  ^  l^2VtV.  In  +  (V'Z  +  V/)  ln(l  -  .  (3-3) 

When  h  ^  A 

E,{h)^ReVaV(,e-^^'^.  (3-4) 

In  this  regime,  where  only  a  negligible  charge  adjustment  is  needed  to  main¬ 
tain  Vs  and  Vb  at  their  asymptotic  values,  the  force  between  the  s  and  b 
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Figure  3-1.  A  spherical  bacterium  of  radius  R  a  distance  h  from  a  wall  surface  (s) 
as  an  idealization  of  a  cell  surface  model.  More  realistic  representa¬ 
tions  of  a  cell  surface  are  shown  in  Figures  3-2  and  3-3. 
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CD  CP  PM  Gl 


Figure  3-2.  Gram-positive  bacterium  surface.  The  lettered  components 
are  explained  in  Reference  [2]. 
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surfaces  is  repulsive  as  long  as  and  Vb  have  the  same  sign.  For  R  = 
lO^A,  €  =  80,  and  characteristic  Vj  =  — 25mV,  =  — 15mVi 

Ec{h)  S  3  •  10“  ergs  {h  >  A).  (3  -  5) 

However,  when  /i  ■C  A  the  repulsion  generally  changes  to  an  attraction.  In 
this  regime 

{-(K  -  Hfln^  +  (V.  +  Hfln2}  (A  >  A).  (3-6) 

As  long  as  V,  ^  14,  in  the  ^  ^  0  limit,  where  a  large  charge  readjustment 

is  needed  to  maintain  the  14,6  surface  potentials,  there  is  an  attractive  force 
between  the  bacterium  and  the  wall.  The  maximum  repulsive  potential  is 
reached  for  Equation  (3-3)  ai  h  =  ho  with 

r ,  y2-\ 

For  the  typical  surface  potentials  which  gave  Equation  (3-5), 

A.  ~  ~  lO-'A  (3-8) 

and 

£^c(max)  =  Ec{ho)  ~  10"^^  ergs.  (3  —  9) 

The  barrier  height 

£^c(nfiax)  kT  ~  5  •  10“^^  ergs  (3  —  10) 

for  T  ~  room  temperature.  From  Equation  (3-5)  the  distance  of  closest 
approach  to  a  wall  with  fixed  surface  potential  14  by  a  spherical  bacterium 
with  fixed  surface  potential  14  and  Brownian  translational  kinetic  energy 
3  kg  T is 
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if  the  only  relevant  forces  are  the  electrostatic  ones  between  a  smooth  wall 
and  a  spherically  symmetric  smooth-walled  inert  bacterium  (cf.  Figure  3-4). 
Such  a  bacterium  would,  therefore,  not  approach  the  wall  closely  enough 
to  feel  the  attractive  potential  of  Equations  (3-6),  and  (3-5)  would  be  an 
adequate  approximation  for  the  electrostatic  part  of  E{h). 

In  the  above  constant  surface-potential  model,  both  the  wall  and  the 
bacterial  surfaces  are  considered  as  conducting  sheets  whose  potentials  are 
fixed  with  respect  to  a  distant  point  in  the  mildly  conducting  fluid  that 
separates  them.  An  alternative  model  makes  the  same  assumption  about  the 
surfaces  but  fixes  the  total  charge  on  each  at  its  asymptotic  value  {h  — ♦  oo) 
rather  than  the  potential.  In  this  case.  Equation  (3-3)  is  changed  to  [4] 

^?(H,  K;  h)  =  E<,(H,  -K;  h).  (3  -  12) 

The  very  short-range  (h  <C  A)  attraction  of  Equation  (3-6)  now  becomes 
a  repulsion,  but  the  long-range  {h  >  A)  repulsion  of  Equation  (3-3)  is  un¬ 
changed. 

Both  models  assume  that  surface  charges  are  free  to  move  along  the  bac¬ 
terial  and  the  wall  surfaces  and  so  cancel  any  tangential  electric  fields:  both 
surfaces  are  described  as  if  they  were  metallic.  But  a  no  less  plausible  model 
might  be  one  in  which  the  surface  charges  do  not  have  such  mobility  during 
the  short  time  of  close  approach  to  the  /tjnjn  of  Equation  (3-11).  Instead,  it 
could  be  assumed  that  the  wall  and  bacterium  surface  charges  remain  fixed  at 
their  asymptotic  positions.  Then  Ec{h)  would  differ  qualitatively  from  that 
of  Equations  (3-3)  or  (3-12)  when  /i  <C  A,  but  it  is  again  given  by  Equation 
(3-4)  when  h  ^  A.^ 

^The  next  term  in  the  expansion  of  Equation  (3-4), 

E2  =  -(H"-fF/)e=^, 
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Figure  3-4.  Electrostatic  interaction  energies  with  constant  potentials  for  different  values  of 
V3andVb  =  -20  mV.  A  -  a. 
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Therefore,  although  the  short  distance  wall-bacterium  Coulomb  inter¬ 
action  is  so  model  dependent  that  quantitatively  reliable  estimates  of  it  are 
questionable,  its  value  ai  h  >  A  and  certainly  for  h  ~  ^min  should  be  un¬ 
ambiguous. 


3.2.2  van  der  Waals  Forces 


The  van  der  Waals  interaction  between  a  bacterium  and  a  wall  can 
become  the  dominant  contributor  to  e{h)  at  distance  h  ^  A  where  Debye 
screening  strongly  suppresses  static  Coulomb  interactions.  The  zero-point 
fluctuations  responsible  for  the  van  der  Waals  interaction  are  mainly  high- 
frequency  ones  (~  10^®  Hz).  The  zero  frequency  screening  of  static  electric 
fields  by  ion  displacement  is  no  longer  relevant  at  these  frequencies. 


The  van  der  Waals  interaction  between  a  pair  of  otherwise  isolated 
molecules  is  always  attractive.  Quite  generally  the  pair  interaction  energy 

£i,  =  Ci,  >0,  (3  -  13) 

as  long  as  the  separation  between  them,  Vij,  satisfies  rij  <  ojq^.  (Otherwise 

Eija  —  Cijc/rJj(jJo)-  Here  is  the  minimum  energy  of  the  excited  states  that 

contribute  to  the  molecular  dipolar  polarizability.  The  total  van  der  Waals 

is  the  sum  of  the  attractive  contribution  from  the  image  force  of  s-charges  with  their 
images  in  the  metallic  6-surface  and  vice-versa.  When  the  surface  charges  are  individually 
fixed  so  that  dielectric  polarization  is  the  only  allowed  charge  movement  in  the  wall  and 
the  bacterium 


where  and  £3  are  the  dielectric  constants  of  the  wall  material,  the  bacterium  surface 

material,  and  the  intervening  liquid,  respectively.  We  note  that  the  image  force  sign 
now  depends  upon  the  differences  between  the  relevant  dielectric  constants  and  may  be 
repulsive  or  attractive. 
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interaction  between  two  condensed  matter  objects  is  often  computed  as  the 
sum  over  all  interacting  pairs  (often  a  quite  incorrect  assumption).  Then  the 
van  der  Waals  attraction  between  two  parallel  walls  separated  by  a  medium 
with  dielectric  constant  €3  would  be 


EM  =  - 


A 

I2vh? 


(3  -  14) 


with 

A  =  nx  n2  Cu  (3  —  15) 


and  with  €3  evaluated  at  some  characteristic  frequency  ~  lOo  and  h  ^ 

The  coefficient  A  is  called  the  “Hamaker  constant”.  However,  it  is  only 
constant  as  long  as  h  ^  When  h  exceeds  A  decreases  as  {cjhijJo). 

For  a  bacterium,  this  decrease  will  begin  when  h  exceeds  several  -lO^A. 


Insofar  as  the  sum  over  pairs  is  a  satisfactory  approximation  for  Ev{h) 
the  van  der  Waals  attraction  between  a  sphere  and  wall  with  the  geometry 
shown  in  Figure  3-1  is 

ATf 

E^h)  =  {h  <  cw„-').  (3  -  16) 

The  Hamaker  constant  for  the  attraction  between  a  bacterium  in  water  and 
a  glass  wall  has  been  estimated  as 

A  ~  5  •  10"^^  ergs.  (3  -  17) 


Summing  Equation  (3-13)  over  all  pairs  can  be  qualitatively  misleading. 
This  is  apparent,  for  example,  if  we  apply  that  procedure  to  the  interaction 
between  bodies  containing  weakly  bound  electrons,  and  especially  to  electron 
conductors  where  oiq  —>■  0.  For  harmonic  oscillators  consisting  of  electrons 
with  mass  m  and  charge  e  bound  with  a  free  oscillation  angular  frequency  ojq 
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when  r  C  cwq  and 


a-i  ~ 


mujQr 


(3  -  19) 


when  r  cwq  ^  so  that  retardation  is  important.  In  both  regimes,  Cij  and  the 
pair-sum  estimate  for  A  diverge  as  u)o  0-  This  would  suggest  an  infinitely 
strong  attraction  between  metallic  plates  (a;o  =  0  for  free  electrons).  This 
incorrect  result  comes  from  the  additivity  assumption.  (It  is  akin  to  using 
an  Einstein  model  instead  of  a  Debye  one  in  considering  the  low-frequency 
excitations  of  a  crystal.)  A  correct  general  way  of  including  all  relevant  effects 
has  been  given  by  Dzyzloshinsii,  Lifshitz,  and  Pitaevskii  (DLP)  [8]. 

There  are  two  families  of  normal  modes  for  electromagnetic  fields  in 
a  homogeneous  medium  with  a  dielectric  constant  e(a;).  One  consists  of 
longitudinal  modes  (e.g.,  plasma  oscillations)  at  those  frequencies  u  =  u^i 
at  which  e{u>i)  =  0.  The  other  is  that  of  the  transverse  modes  for  which 
P  =  For  the  inhomogeneous  geometry  of  Figure  3-5,  the  nor¬ 

mal  mode  frequencies  depend  upon  h  and  the  three  dielectric  constants 
€i(a;),c(a;),  €3(0;).  Each  normal  mode  with  frequency  w„(/i)  has  a  zero-point 


energy 


En  = 


hu}n{h) 


Then  the  force  between  the  walls  1  and  2  of  Figure  3-5  is 

^  =  -dh\'~T- 

and  the  /i-dependent  part  of  the  zero-point  energy  is 

Ev{h)  =  23  ^K(h)  -  u;„(oo)]. 


(3  -  20) 


(3-21) 


(3  -  22) 


DLP  find  that  the  longitudinal  mode  contribution  to  Equation  (3-22)  is  dom¬ 
inant  at  sufficiently  small  h  and  the  Hamaker  constant 


h  <  cjoj 


(3  -  23) 
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^  =  r  .  (3  -  24) 

Jo  [cl(^a)  +  C3(ta)J  [C2(m)  —  €3(ia}J 
(The  ratios  in  the  brackets  are  just  those  that  enter  into  the  image  forces 
from  each  surface.)  The  analytic  continuation  to  imaginary  frequency  in  the 
€(ia)  of  Equation  (3-24)  is  easily  accomplished  and  related  to  measured  e(a;) 
through  the  Kramers- Kronig  relation.  This  states  that  any  linear  causal 
system  (say,  of  electrons)  gives  the  same  response  to  an  electric  field  as 
that  from  a  sum  of  infinitesimally  damped  harmonic  oscillators  with  various 
frequencies  (w') 

rjW.  (3_25, 

'  '  Jo  (u?')^  — 

The  number  density  of  such  oscillators  with  free  natural  frequencies  between 
(jJi  and  (j}dwi, p{u}Wi)dwi  satisfies  the  sum  rule 


(3  -  25) 


f  p{u)')du}'  = 
Jo 


4jrne^ 


(3  -  26) 


Here  n  is  the  spatial  number  density  of  electrons.  (We  assume  that  virtual 
electron  transitions  dominate  the  van  der  Waals  attraction  between  bacteria 
and  walls.)  From  Equation  (3-25)  with  u)  approaching  the  real  axis  from 


above 


-Im  c(u>)  =  p{ijt}) 


(3  -  27) 


/  •  X  ,  2  r°°  (jj  Im  c(a;) 

e(ia)  =  14--/  — y-du?. 

“K  Jo  (XT 


(3  -  28) 


The  e(m)  of  Equations  (3-28)  and  (3-24)  are  positive,  monotonically  decreas¬ 
ing  functions  of  a  that  can  be  quantitatively  calculated  whenever  sufficient 
measurements  of  absorptive  properties  {Im  c(w))  of  materials  1,  2,  and  3  are 
known.  ^ 

^The  needed  correction  to  the  wq  — ^  0  divergence  in  A  when  Equations  (3-18)  and 
(3-19)  are  inserted  into  Equation  (3-15)  is  now  easily  seen.  Consider  materials  1  and  2  to 
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Several  general  properties  of  the  Hamaker  constant  follow  from  Equa¬ 
tions  (3-23)  and  Equation  (3-24). 

(a)  >  0  (van  der  Waals  force  is  attractive)  if  media  1  and  2  have 
the  same  dielectric  constants. 

(b)  A  >  0  (i.e.,  a  bacterium-wall  van  der  Waals  “attraction”  becomes 
a  repulsion)  if  appropriately  averaged  e{iu})  satisfy 

Cl  >  C3  >  C2  a) 


or 

Cl  <  €3  <  €3.  b)  (3  —  29) 

[This  is  the  case,  for  example,  if  1  is  quartz,  2  is  decane,  and  3  is 
albumen.] 

(c)  When  Ci  =  cj,  A  is  unchanged  if  medium  3  is  interchanged  with 
both  medium  1  and  medium  2. 

(d)  When  (ci,2  —  1)  >  (ca  —  1)  the  Hamaker  A  is  diminished  by  the 
presence  of  medium  3  relative  to  its  value  when  C3  =  1. 


A  variety  of  computed  values  of  the  Hamaker  A  (neglecting  retardation 
corrections)  based  upon  measured  values  of  Im  e(a;)  are  given  in  Tables  3.1 
and  3.2,  together  with  some  A  computed  with  Equation  (3-15)  [3].  The 


consist  of  electron  harmonic  oscillators  with  frequency  uq  so  that 


fi,2  =  1  ~ 


and  assume  fa  =  1  (or  equivalently,  take  Ci,2  =  1  and  for  (3  the  RHS  of  fi,2  above).  Then 
from  Equations  (3-23)  and  (3-24) 


In  the  met^LIhc  limit,  wq  0,  and  A  0.13u;p .  Quite  generally,  summing  over  “isolated” 
molecular  pairs  to  approximate  A  is  valid  only  if  wo  ^Wp. 
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Table  3.1.  Computed  values  for  the  Hamaker  constant  A  for  identical  walls 
separated  by  water.  Equations  (3-23)  and  (3-24)  are  used  with  measured  Im  €{u) 
and  Equation  (3-28)  except  for  values  marked  *  where  a  version  of  the  summed 
pair  approximation  is  made  [3]. 


Walls  (separated  by  H2O) 

A  (10  ergs) 

Cl, 2  =1 

26 

decane 

3* 

benzene 

2 

borine  albumen 

7 

dipalmytol 

5* 

quartz 

19* 

mica 

15* 

CS2 

10 

Table  3.2.  Same  as  Table  3.1  except  the  medium  between  the  walls  has  Ca  =  1 
[3]. 


Walls  (separated  by  62  =  1) 

A  (10  ergs) 

water 

26 

C2S 

82 

decane 

54* 

quartz 

100* 

mica 

15* 
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Table  3.3.  Retardation  correction  factors  for  Hamaker  A  [3]. 


1  —  3  —  2 

0 

decane  -  vacuum  -  decane 

0.3 

decane  -  water  -  decane 

0.2 

dipalmitoyl  -  water  -  dipalmitoyl 

0.3 

mica  -  water  -  mica 

0.2 

(ci)  =  1  -  water  -  (62)  =  1 

0.4 

borine  albumin  -  water  -  borine  albumin 

0.3 

quartz  -  ablumin  -  decane 

0.08 

intervening  medium  is  water  in  Table  3.1  and  unpolarizable  in  Table  3.2. 
Medium  1  is  always  the  same  medium  2. 

When  all  three  media  are  different,  the  computed  A  may  be  qualitatively 
altered.  For  example,  a  quartz  and  a  decane  wall  separated  by  albumin  give 
a  repulsive  A  =  —1.3  •  10“^^  ergs. 

When  k  >  cu>~^  retardation  effects  must  be  included.  These  reduce  the 
effective  by  a  factor  of  order  c/hH).  Calculated  ratios  for  ^  =  i4(500A)/>l(10A) 
are  given  for  various  1-3-2  wall-liquid-wall  combinations  in  Table  3.3.  Because 
separations  of  interest  are  often  several  lO^A  this  correction  factor  is  quan¬ 
titatively  important. 

Over  2000  calculations  have  been  made  [3]  for  A  when  walls  1  and  2 
are  those  of  a  cell  and  medium  3  is  water  (Figure  3-6).  Calculations  were 
made  for  various  assumed  surface  cost  thicknesses  and  compositions.  All  the 
computed  A  values  (at  h  =  50  k)  are  between  1  •  10"^^  ergs  and  4*  10~*^  ergs, 
with  A(200A)/A(50A)  0.3  1.0.  The  assumed  glass- water-cell  Hamaker 

constant,  A  ~  5  •  10"^^  ergs,  near  the  geometric  mean  for  cell-water-cell  and 
quartz-water-quartz,  would  seem  plausible. 
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Interior  of  Cell 


Figure  3-6.  Model  for  a  cell  surface  used  in  van  der  Waals  force  calucations  [3]. 
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3.2.3  Wall-Cell  Attraction  at  Large  Separations 


From  Equations  (3-15),  (3-16),  and  (3-17) 

E{k)  =  E,(h)  +  E,(h)  =  3  •  10-"  (^5^)  e"*'*  -  8  ■  lO"''  (jj  ergs 

(3  -  30) 

for  A  ■C  A  <  Without  the  van  der  Waals  attraction,  the  distance  of 

closest  approach  of  the  rigid  spherical  bacterium  and  a  wall  was  the  ~ 

7A  of  Equation  (3-11)  where  Ec{h  min)  ~  ksT.  At  that  distance,  the  van  der 
Waals  attraction 

•B.(Amin)  ~  8  ’  10"“^  ~  2  '  lO’"  (j)  ksT.  (3  -  31) 

Then  E{h)  could  be  attractive  and  exceed  ksT  only  if 

A  <  2  •  lO-^i?  ~  2  •  lO^A,  (3  -  32) 

corresponding  to  an  electrolyte  concentration  in  excess  of  about  10“^^  moles/litre. 

The  van  der  Waals  attraction  should  not  be  effective  in  holding  the  mod¬ 
eled  bacterium  near  a  wall  unless  A  is  very  much  less  than  that  in  Equation 
(3-32).  The  minimum  time  it  would  take  for  Brownian  motion  to  remove  a  re¬ 
versibly  attached  bacterium,  modeled  as  described  above,  from  an  approach 
distance  to  a  wall  of  is 

fSiiniall 

•‘bt  J  (3-33) 

with  the  bacterium  diffusion  constant  (momentum  x  velocity^ /Stokes  drag) 

D  ~  10“®cm^s“^.  (3  —  34) 

For  van  der  Waals  attachment  to  last  for  t  >  lO^s,  ■E’(^niin)  — 

Therefore,  even  in  the  absence  of  any  electrostatic  repulsion  between  the 
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bacterium  and  a  (silicate)  wall,  an  ~  5-10“^'*  ergs  Hamaker  constant  would 
give  a  significant  attachment  only  if  h  <  lO^A.  When  electrostatic  repulsion 
is  included,  an  assumed  rigid  spherical  Achromobacter  R8  bacterium  {R  ~ 
OAfi)  could  attach  with  van  der  Waals  forces  alone  for  more  than  10^  s  only  if 
A  <  20A  <  80A).  However,  there  is  evidence  that  Achromobacter  R8 

bacteria  can  accomphsh  strong  attachment  even  when  A  >  20 A  [9].  Fig  ure 
3-7  shows  how  the  measured  attachment  after  an  hour  (No)  varies  with  A 
[10].  The  number  of  attached  bacteria  is  well  represented  by 


Na  =  No  In  A  <  Ao  a) 

=  0,  A  >  Ao  b) 


(3  -  35) 


with  Ao  ~  250A.  With  this  Debye  length,  the  van  der  Waals  attraction  at 
the  minimum  distance  of  approach  is  only  about  ksT,  and  it  is  difficult  to 
see  how  it  would  play  a  major  role  in  near-threshold  attachment  (A  ~  Aq)  if 
the  bacterium  is  adequately  described  by  the  model.^ 

Therefore  modelling  a  bacterium  as  an  inert  spherically  symmetric  charged 
shell  interacting  with  a  similarly  charged  smooth  conducting  wall  seems  much 
too  simple  to  describe  bacterium-wall  adhesion  or  even  how  the  bacterium 
attaches  temporarily  (and  reversibly)  for  a  time  long  enough  (>  lO^s?)  for  a 
polymer  bridge  to  grow  between  the  bacterium  and  the  wall. 


One  obvious  oversimplification  is  the  assumption  of  exact  spherical  sym¬ 
metry  for  the  bacterium.  In  an  electrolyte,  the  electrostatic  interaction  be¬ 
tween  a  bacterium  and  a  wall  falls  off  so  rapidly  with  separation  h  in  the 

attachment  is  accomplished  by  one  or  more  surface  molecules  reaching  out  beyond 
the  nominal  spherical  surface  to  attach  to  the  wall,  Equation  (3-35)  would  imply  that  the 
probability  for  such  an  extension  to  a  distance  between  /  and  I  +  dl  is  P(l)dl  =  kdl/l  if 
I  <lo  and  zero  if  /  >  /<,,  with  ~  lO^A.  This  is  very  different  from  the 


logP(/)a- 


{i-ionf 


/-distribution  of  a  randomly  oriented  chain  of  length  /«  and  link  length  A. 
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Figure  3-7.  Effect  of  electrolyte  concentration  on  A  an  Achromobacter  R8  attachment 

to  a  wall.  The  electrolytes  are  Na+Cl"  and  Mg‘*"*‘S04“  “  (by  K.  Marshall  [10]). 
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A  A  regime  that  only  the  parts  of  a  bacterium’s  surface  charge  nearest  the 
wall  are  relevant.  Similarly,  because  of  the  fall-off  in  Equation  (3-13), 
only  the  part  of  the  bacterium  surface  nearest  the  wall  is  relevant  in  Equation 
(3-16).  Therefore,  in  Equation  (3-30),  the  radius  R  in  the  first  term  on  the 
RHS  is  the  radius  of  curvature  of  the  bacterium’s  surface  charge  layer  and  in 
the  second  is  the  radius  of  curvature  of  its  polarizable  surface  matter,  both 
to  be  evaluated  at  the  position  on  the  surface  nearest  the  wall.  There  seems 
to  be  no  compelling  reason  for  them  to  be  everywhere  equal  to  an  average 
radius  or  to  each  other  in  a  living  bacterium.  Then,  if  some  areas  of  the 
bacterial  surface  are  flatter  in  their  distribution  of  polarization  matter  with 
a  local  radius  of  curvature  >  R,  the  van  der  Waals  attraction  would  be 
increased  by  Rc/R  when  the  bacterium  orientation  brings  that  area  closest 
to  the  wall.  A  similar  increase  could  come  from  microscopic  wall  irregular¬ 
ities  as  long  as  those  parts  of  the  wall  nearest  the  bacterium  do  not  carry 
wall  surface  charge,  ff  14,6  ^  —2QmV  and  A  ~  lO^A,  the  average  separation 
between  surface  charges  on  both  objects 

Surface  irregularities  might  then  bring  uncharged  parts  of  the  wall  much 
nearer  the  bacterium,  where  the  combination  of  increased  attractive  van  der 
Waals  force  and  an  increased  attractive  image  force  (cf.  second  equation  in 
footnote  1)  might  be  enough  to  allow  reversible  binding  during  which  a  more 
permanent  colloidal  bridge  could  cement  cell-wall  attachment. 

If  the  general  shape  of  a  bacterium  departs  greatly  from  a  spherical  one, 
the  wall-bacterium  van  der  Waals  attraction  can  become  much  larger.  Thus, 
if  the  sphere  of  radius  R  is  replaced  by  a  cylinder  of  radius  R  and  length 
L,  then  E^{h)  depends  sensitively  upon  the  orientation  of  the  cylinder.  It 
is,  or  course,  largest  when  the  cylinder  axis  is  in  a  plane  parallel  to  the  wall 
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surface.  In  this  case,  the  attractive  van  der  Waals  potential  of  Equation 
(3-16)  is  increased  by  a  factor  {Lf  Rhy^^.  This  could  be  sufficient  to  make 
Ey(h)  much  more  than  kaT  even  at  A  ~  lO^A.  If  the  bacterium  could  offer 
a  flat  surface  to  the  wall,  the  increase  in  E„(/i)  would  be  even  greater. 

The  very  complex  and  dynamic  structure  of  a  hving  cell  may  allow 
probes  and  protuberances  almost  to  touch  a  wall  even  when  the  average 
distance  of  closest  approach  ^min  ~  7A  ~  lO^A.  If  the  effective  radius  of 
an  uncharged  “probe”  head  exceeds  30  times  its  approach  distance,  then 
van  der  Waals  attraction  can  attach  it  to  a  wall  with  a  binding  energy  > 
lOksT.  A  single  chemical  bond  or  ion-image  binding  could  also  give  a  similar 
attachment  energy. 

Finally,  natural  walls  are  very  far  from  smooth  on  the  sub-micron  scale. 
Their  representation  as  smooth  uniform  surfaces  ignores  special  properties  of 
very  small  patches  that  could  increase  the  strength  of  attraction  to  certain 
cell  surface  areas  as  well  as  pervasive  irregularities. 

Groundwater  slowly  moving  through  rock  pores  and  capillaries  would  be 
expected  to  have  a  large  electrolyte  concentration:  if  very  old  it  may  exceed 
10"^  m  which  would  give  A  <  ISA.  In  such  strong  electrolytes  Equation 
(3-31)  gives  >  lAk^T  at  /i  ~  7A  <  lO^A  so  that  semi- permanent,  but 
reversible,  van  der  Waals  attachment  to  walls  may  be  expected  even  without 
the  help  of  ultimate  polymer  bridging.  In  the  strong  electrolyte  regime,  it 
does  not  appear  necessary  to  understand  the  complicated  active  nature  of 
a  cell  surface  to  account  for  expected  attachment.  However,  penetration 
of  a  spherically  symmetric  bacterium  actually  to  touch  a  smooth  uniformly 
charged  wall  surface  would  still  require  overcoming  a  barrier  of  order  lO^ksT. 
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The  attractive  force  for  the  van  der  Waals  interaction  of  a  spherical 
bacterium  with  a  smooth  glass  wall  (F)  is  about  E{h)lh: 


lO^A 


dynes. 


(3-37) 


This  force  cannot  be  compared  directly  to  that  from  possible  van  der  Waals 
attachment  because  the  detachment  left  “polymer  footprints”  on  the  attach¬ 
ing  wall  surface:  the  polymer  bridges  that  formed  in  attachment  and  were 
broken  in  detachment  may  have  contributed  the  major  attachment  force. 
The  force  from  a  single  bacterial  molecule  held  to  a  wall  by  covalent,  ionic, 
or  image  force  binding  would  contribute  a  detachment  threshold  force 


2  •  10  ^  dynes 


(3  -  38) 


about  equal  to  the  largest  measured  one  [10]. 


The  importance  of  van  der  Waals  forces  in  attachment  might  be  explored 
by  experiments  in  which  wall  materials  are  varied.  If  a  wall  is  made  of 
material  that  satisfies  Equation  (3-29),  the  van  der  Wcials  interaction  is  no 
longer  attractive.  [If  a  glass  or  quartz  wall  could  be  replaced  by,  say,  a  copper 
one  with  a  sufficiently  thin  (<C  10^ A)  glass  or  quartz  surface  layer,  the  van 
der  Waals  part  of  the  cell-wall  interaction  with  a  water  intermediary  might 
be  greatly  diminished  or  reversed.] 


3.3  Assessment  of  the  Energetics  and  Kinetics  of  Ad¬ 
hesion 


In  the  context  of  bioremediation,  the  major  motivation  for  developing 
quantitative  measures  of  bacterial  adhesion  is  to  be  able  to  predict  rates  of 
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transport  in  porous  media.  Unfortunately,  as  described  in  the  previous  Sub¬ 
section  (3.2),  the  theoretical  basis  for  understanding  adhesion  of  bacteria  to 
walls  is  limited  to  the  Derjaguin-Landau-Verwey-Overbeek  (DLVO)  model 
developed  to  describe  adsorption  of  collioids.  While  this  model  allows  us  to 
understand  the  important  physical  factors  in  adhesion  qualitatively,  it  cannot 
be  used  for  quantitative  prediction  of  the  behavior  of  any  real  bacterium.  Fur¬ 
thermore,  with  a  few  exceptions,  experimental  studies  of  bacterial  adhesion 
available  in  the  literature  provide  only  qualitative  information.  There  are 
very  few  qualitative  measurements  of  energetics  or  kinetics  of  adhesion  that 
could  be  used  to  evaluate  rates  of  mass  transport,  even  for  specific  systems. 
Studies  that  would  allow  general  quantitative  predictions  relating  adhesion 
to  cell-wall  properties  are  nonexistent.  Thus,  the  analysis  we  present  below 
is  made  for  a  generic  case  to  illustrate  the  nature  of  the  information  needed 
and  the  ranges  of  parameters  that  are  physically  reasonable. 

The  DLVO  model  of  cell-wall  adhesion  can  be  used  to  illustrate  the  gen¬ 
eral  trends  of  bacterial  adhesion.  The  interaction  strength  can  be  computed 
as  a  function  of  distance  between  the  bacterium  and  the  wall  using  Equations 
(3-2),  (3-3),  and  (3-16)  with  representative  values  for  the  effective  potentials 
and  the  Hamaker  constant  [3],  [11],  [12].  The  results  are  shown  in  Figure  3-8 
for  a  range  of  solution  ionic  strengths  comparable  to  sea  water  (2  x  10“^  M), 
old  groundwater  (2  x  10"^  M),  and  new  groundwater  (2  •  10"^  M).  It  is  imme¬ 
diately  apparent  that  the  nature  of  bacterial  adhesion  changes  qualitatively 
over  the  range  of  ionic  strengths  relevant  to  subsurface  remediation.  At  high 
ionic  strengths,  the  adhesive  interaction  is  strongly  attractive,  resulting  in 
binding  energies  much  larger  than  the  thermal  energy,  and  thus  irreversible 
adhesion  [13]  [14].  At  low  ionic  strengths,  the  repulsive  electrostatic  interac¬ 
tion  is  so  strong  that  adhesion  is  impossible  if  only  the  physical  (as  opposed  to 
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Figure  3-8.  Interaction  energy  as  a  function  of  distance  from  the  wall  calculated  using  the 
DLVO  model  described  in  Section  3.2.  The  curves  show  the  charactenstic 
secondary  minimum  in  potential  energy,  which  is  separated  by  a  strong  potential 
energy  barrier  from  the  primary  minimum  (not  shown)  at  less  than  lOA  fron^e 
surface.  Parameters  used  in  the  calculation  were  a  bacterial  racKus  of  l|xm,  30mV 
and  20mV  for  the  effective  potentials  on  the  bacterium  and  wall  respectivdy. 
Values  of  the  distance-dependent  Hamaker  constants.  A,  used  for  calculating 
cell-wall  interactions  are  shown  below. 
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biological  chemical)  mechanism  of  adhesion  is  in  effect.  The  tremendous  vari¬ 
ability  in  the  strength  of  interaction  with  ionic  strength  has  two  important 
consequences.  The  first  is  that  bacterial  adsorption  may  change  dramatically 
if  the  environment  changes  as  the  bacteria  move  through  the  subsurface.  The 
second  is  that  model  predictions  of  the  interaction  are  so  sensitive  to  the  pa¬ 
rameters  chosen  that  experimental  measurements  are  essential  to  accurate 
prediction. 


Once  the  interaction  strength  is  known,  then  the  kinetics  of  adhesion 
can  be  determined  using  standard  approaches.  The  motion  of  a  bacterium 
(in  the  absence  of  biological  motor  motion)  is  described  by  the  Langevin 
equation  [15] 

dv  ,  ,  ,  . 

^  ^waU  -  +  F(t)  (3  -  39) 

where  i^waU  is  the  force  due  to  the  interaction  with  the  wall,  which  is  opposed 
by  the  viscosity  of  the  water  rj  acting  on  the  motion  of  the  spherical  bacterium 
of  radius  rtacj  and  F{t)  is  the  stochastic  force  on  the  bacterium  due  to  the 
random  colhsions  of  surrounding  water  molecules.  Far  from  the  wall,  the 
bacterium’s  motion  will  be  Brownian,  with  a  diffusion  coefficient 

kT 

D  =  - - =  2.2  X  10"®  cm^s"^  (3  -  40) 

07r7^7*l)ac 

for  a  radius  of  Ifim  and  the  viscosity  of  water  t]  =  0.01  gfcm^.  In  a  region  of 
constant  force,  the  bacterium  will  be  accelerated  to  a  steady  state  velocity. 


F 

Ott  J/rbac 


(3-41) 


The  time  constant  to  reach  steady  state  is  ti/2  =  0.1  ms,  which  allows  a 
random  Brownian  displacement  of  approximately  100  A  in  the  time  required 
to  reach  1/2  of  the  steady  state  velocity. 
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A  qualitative  analysis  of  the  rates  of  adsorption  based  on  these  concepts 
is  suflScient  for  the  purpose  of  illustration  here.  The  forces  near  the  wall 
for  the  interactions  in  intermediate  ionic  strength  solution  (solid  curve  in 
Figure  3-8)  are  shown  in  Figure  3-9.  The  bacterium  experiences  repulsive 
forces  on  the  order  of  10“®  dynes  near  the  wail,  as  shown  in  Figure  3-9a. 
However,  during  an  approach  from  solution,  it  first  experiences  attractive 
forces  pulling  it  into  the  secondary  minimum  of  potential  energy.  These  forces 
are  two  orders  of  magnitude  smaller  than  the  repulsive  forces  near  the  wall, 
but  they  are  still  sufficient  to  accelerate  the  bacterium  to  speeds  of  lOfim/s 
as  it  approaches  the  secondary  minirrium.  A  quantitative  analysis  of  the 
kinetics  of  approach  to  the  wall  using  Equation  (3-39)  will  yield  a  macroscopic 
rate  parameter  governing  adsorption.  This  is  most  typically  described  as  a 
sticking  probability  or  sticking  coefficient,  which  is  just  the  probability  that 
a  particle  approaching  the  walls  adsorbs  rather  than  returning  to  the  bulk 
solution.  From  the  estimates  of  the  speed  of  approach,  it  is  clear  that  in 
the  moderate  ionic  strength  solution,  the  stiking  probability  will  be  close  to 
one.  However,  with  decreasing  ionic  strength,  the  probability  will  decrease 
rapidly,  until  in  the  case  of  low  ionic  strength,  shown  by  the  long-dashed 
curve  of  Figure  3-8,  it  will  be  essentially  zero. 

Once  a  sticking  probability  is  determined,  the  rate  at  which  bacteria  will 
adhere  to  walls  can  be  calculated  using  standard  methods  of  mass  transport, 
without  having  to  incorporate  the  potential  energy  near  the  wall  into  the 
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Figure  3-9.  Force  exerted  by  the  wall  on  a  bacterium  as  a  function  of  distance  from  the  wall. 

The  values  are  calculated  for  the  interaction  potential  of  moderate  ionic  strength 
(solid  curve)  of  Figure  3-8.  The  steady  state  speeds  corresponding  to  the  forces 
were  calculated  using  Equation  (3-4  la.),  a)  Strong  repulsive  forces  near  the  wall  are 
dominated  by  the  strong  mawmum  in  the  potential  energy  at  approximately  30A. 
b)  At  larger  distances  weakly  attractive  forces  accelerate  approaching  bacteria  into 
the  secondary  minimum  of  the  potential  energy  curve. 
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analysis.'^  The  problem  of  interest  is  the  case  where  a  clean  wall  is  suddenly 
exposed  to  a  solution  of  bacteria  of  concentration  Cq.  Mass  transport  is 
governed  by  the  equation 

>(!,()« (3-42) 

where  j  is  the  flux  to  the  wall,  c  is  the  concentration,  and  x  is  the  dis¬ 
tance  from  the  wall,  which  must  be  combined  with  the  requirement  of  mass 
conservation, 

dc{x,t)  _  djix,t)  _  . 

dt  dx  '  ^  ^ 

The  effect  of  the  attractive  potential  due  to  the  wall  is  accounted  for  by  the 

boundary  conditions,  which  are  expressed  in  terms  of  the  sticking  probability 

p,  as 

c(0,t)  =  (1  —  p)co,  and  c(oo,t)  =  Cq.  (3  —  44) 


This  is  a  standard  problem  in  adsorption,  and  yields  the  solution  for  the 
flux  to  the  wall  [16]: 

This  has  the  physical  meaning  of  a  drop  in  concentration  from  the  bulk 

solution  value  of  c(oo)  =  Cq  to  c  =  (1  “p)co  over  a  diffusion  layer  of  thickness 

the  potential  energy  is  included  explicitly,  the  one-dimensional  flux  of  bacteria 
toward  the  wall  must  be  calculated  using  the  standard  expression 

.  _  —Dc  d/i 

^  ~  'W  d^’ 

where  c  is  the  concentration  of  bacteria  and  fi  is  the  chemical  potential,  which  includes  a 
term  due  to  the  potential  energy  V(x)  near  the  wall, 

//  «  fcT  In  c-h  V{x), 


yielding 


dc  -Dc 
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\/irDt.  Using  the  value  of  the  diffusion  coefficient  from  above,  and  a  solution 
concentration  of  10®/cm^,  we  find  the  rate  of  absorption  of  bacteria  onto  the 
wall  in  any  time  interval  to  be  approximately 

AN{t2  —  h)  =  p(5  X  lO^/sl/2  —  ca:?){ty^  —  t^^).  (3  —  46) 

The  fraction  /  of  the  bacteria  lost  from  solution  in  a  tube  of  radius  r  then 
becomes 

/  =  -(5  X  10"®cm/s^/^)(t2^^  -  (3  -  47) 

T 

This  rate  of  loss  would  continue  until  the  density  of  adsorbed  bacteria  became 
large  enough  to  block  adsorption  of  additional  bacteria,  and  then  would  fall 
to  match  the  rate  of  desorption  when  the  equilibrium  surface  concentration 
was  reached.  The  effects  of  rates  of  loss  of  this  magnitude  on  transport 
through  porous  media  are  illustrated  in  Subsection  3-4  for  a  case  equivalent 
to  a  sticking  probability  of  p  =  10"^. 

From  the  discussion  presented  above,  it  should  be  clear  that  a  realistic 
assessment  of  the  rate  of  adsorption  of  bacteria  onto  walls  (and  thus  the 
rate  of  transport  through  porous  media)  cannot  be  made  in  the  absence  of 
experimental  measurements.  Somewhat  surprisingly,  there  is  very  little  data 
available  in  the  literature  to  provide  the  needed  information.  For  specific 
systems,  standard  batch  measurements  in  which  the  number  of  bacteria  lost 
from  solution  is  monitored  as  a  function  of  time  would  be  useful  if  care¬ 
fully  designed.  However,  the  results  of  such  measurements  are  difficult  to 
extrapolate  to  variable  conditions.  A  better  alternative  would  be  to  perform 
experiments  that  measure  the  forces  of  adhesion  directly.  Such  measurements 
would  yield  information  that  can  be  directly  applied  to  predictions  of  kinet¬ 
ics  and  that  would  allow  the  nature  of  the  mechanisms  to  be  determined  so 
that  extrapolations  to  different  conditions  could  be  made  with  confidence.  A 
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major  question  that  should  be  addressed  by  such  studies  is  under  what  con¬ 
ditions  physical  interactions  dominate  and  under  what  conditions  biological 
interactions  dominate. 

The  apphcation  of  direct  force  measurements  to  the  study  of  adhesion  is 
thoroughly  understood  [17].  For  the  study  of  bacterial  adhesion,  at  least  two 
types  of  measurement  techniques  could  be  applied  to  quantify  the  force  on  the 
bacterium  as  a  function  of  distance  from  the  surface.  One  technique  would  be 
to  combine  the  interferometric  microscope  measurements  that  are  already  in 
use  to  study  bacterial  adhesion  [18]  with  laser  trapping  techniques  that  are  of 
increasing  utility  in  microbiology  [19].  The  forces  that  can  be  generated  on  a 
bacterium  with  a  focused  laser  beam  can  be  as  large  as  10“®  dyne  at  powers 
that  do  not  damage  the  bacterium  [20].  This  is  sufficient  to  probe  the  long- 
range  forces  of  chemical  and  physical  interaction  as  illustrated  in  Figure  3-9. 
However,  it  will  not  allow  direct  measurement  of  the  short-range  forces  nor 
of  the  forces  due  to  biological  adhesion  discussed  in  Subsection  3.2.  The  laser 
trap  could  be  used  to  place  a  bacterium  near  the  surface,  and  its  motion  could 
be  followed  dynamically  following  release  to  measure  rates  of  adsorption  [45]. 
Alternatively,  the  trapped  bacterium  could  be  moved  to  different  distances 
from  the  surface  by  moving  the  beam,  and  the  force  of  interaction  determined 
by  varying  the  laser  intensity  until  the  trapping  force  was  insufficient  to  match 
the  force  due  to  the  surface.  Finally,  the  bacterium  could  be  held  at  a  fixed 
position,  and  the  time  constant  for  biological  adhesion  to  occur  and  pull  the 
bacterium  out  of  the  trap  could  be  determined.  This  approach  has  the  benefit 
that  it  is  non-intrusive  for  the  bacterium.  It  has  the  limitations  of  requiring  a 
transparent  and  smooth  substrate  (for  the  interference  microscope),  so  that 
the  role  of  substrate  heterogeneity  could  not  be  probed  readily. 

Another  measurement  technique  that  could  be  used  is  Atomic  Force 
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Microscopy  (AFM)  [21]  [22].  In  this  technique,  the  force  between  an  object 
and  a  surface  is  measured  by  placing  the  object  on  a  very  soft  cantilever 
beam  and  measuring  the  deflection  of  the  cantilever  as  it  is  moved  relative 
to  the  surface.  Forces  as  small  as  10”^  dynes  are  readily  detected  in  this  way 
so  that  short-range  physical  and  biological  interactions  could  be  measured 
directly  using  AFM.  This  technique  would  require  devising  methods  of  at¬ 
taching  a  bacterium  to  the  cantilever  without  perturbing  its  function.  Once 
accomplished,  the  technique  would  allow  adhesion  properties  to  be  measured 
under  any  conditions,  and  in  particular,  the  influence  of  surface  heterogene¬ 
ity  could  be  probed  directly  by  moving  the  probe  to  different  locations  above 
a  rough  surface. 

3.4  Flow  Induced  Force  on  Wall- Attached  Bacteria 


In  the  regime  of  interest,  there  is  a  very  slow  flow  of  liquid  through  very 
narrow  capillaries  and  cracks  to  a  depth  far  below  the  earth  surface.  This 
flow  is  assumed  to  be  driven  by  the  weight  of  the  fluid.  The  average  fluid 
speed  in  a  capillary  with  constant  radius  is 


V  = 


ga^d 


(3  -  48) 


where  rj  =  kinematic  viscosity  of  the  fluid 
a  =  capillary  radius 

d  =  the  change  in  depth  from  top  to  bottom  of  the  capillary 
I  =  the  total  capillary  length 
g  =  the  acceleration  of  gravity. 

For  a  =  10/i,  j/  =  5 •  10“^  (water  at  50°),  and  d/l  =  10~^, v  =  3- 10“^  cm  s“^. 


41 


The  Reynolds  number  of  the  flow 


Re  =  —  <  5  •  10-^(ai3) 
V 


(3  -  49) 


is  far  below  unity  if  a  -C  so  that  the  flow  would,  in  such  cases,  be 

expected  to  be  smooth  (Pouiselle).  At  the  capillary  wall  (Figure  3-10)  the 
fluid  velocity  is  zero  and  its  average  value  near  an  attached  bacterium  of 


radius  R  is 


R  _  gaRd 

Vb - v  = 

a  oTfl 


(3  -  50) 


The  Stokes  drag  on  the  bacterium  from  this  flow  is,  very  approximately, 
^  ^  ddR  /rt  p- 1  \ 


Fb  ~  dirpiTjRvb  ~  pig- 


(3  -  51) 


with  Pi  the  liquid  density.  This  (tangential)  drag  force  is  independent  of  rj. 


With 


Pi  =  ^(HaO)  =  Ig  cm* 


the  tangential  drag  is  then,  roughly. 


h  ~  W,--, 


(3  -  52) 


(3  -  53) 


with  Wb  the  weight  of  the  attached  bacterium.  For  R  ~  1/x,  a  ~  lOR,  pi 
Ig  cm“^,  and  /  ~  lOd  Equation  (3-51)  gives 


^6  ~  10  ^  dynes. 


(3  -  54) 


This  is  a  purely  tangential  drag.  A  lift-off  force  enters  only  when  higher 
orders  in  Vb  are  included  in  calculating  the  force  of  the  flowing  viscous  fluid 
on  the  bacterium.  (If  Re  1,  with  turbulent  boundary  layers,  it  is  expected 
to  be  of  order  Fd.)  For  very  small  Re,  as  is  the  case  here,  any  small  lift 
component  should  not  exceed  the  Oseen  correction  to  the  Stokes  drag: 


^^(lift)  <  Fb- 


PiR^vl  ~  10 


(3  -  55) 
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a  force  which  is  entirely  negligible.  However,  if  a  bacterium  is  kept  from 
sliding  along  a  surface  in  response  to  the  tangential  Fb  of  Equations  (3-51) 
and  (3-52)  by  its  interaction  with  the  surface,  the  flow  induced  depinning 
force  might  approach  these  values.  This  can  occur  in  several  ways: 


(1)  Geometrical  path  bifurcations.  These  can  stop  the  bacterium’s  “glide” 
in  response  to  i>  at  intersection  (stagnation)  points  like  those  at  point 
A  in  Figure  3-11.  The  squeeze  on  the  bacterium  from  the  two-sided 
flow  could  give  a  detachment  force 

Fd  ~  -Fb,  (3  -  56) 

fc 

with  Tc  the  radius  of  curvature  of  the  attaching  surface  at  the  point  A. 

(2)  Tangential  variations  in  surface  properties.  The  attachment  binding 
energy  of  a  gliding  attached  bacterium  (Eb)  can  vary  with  distance  (s) 
along  the  surface.  This  will  slow  the  tangential  motion  of  an  attached 
bacterium  and  cause  a  detaching  torque  from  the  resulting  Stokes  drag. 
If  dEbldg  exceeds  that  Stokes  drag,  the  bacterium  would  not  move 
along  its  attachment  surface  in  response  to  the  flow.  The  effective 
detachment  force  is  at  most  that  of  Equation  (3-55). 

(3)  Bacterium  attachment  to  the  capillary  surface  by  polymers  that  an¬ 
chor  to  specific  sites  on  the  surface.  These  bridge  polymers  keep  the 
bacterium  from  moving  in  response  to  the  Stokes  drag  unless  that  drag 
exceeds  the  yield  strength  of  the  bridge.  However,  the  calculated  at¬ 
tachment  forces  of  Equations  (3-37)  and  (3-38)  greatly  exceed  the  max¬ 
imum  flow  induced  lift. 
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3.5  Conclusions  and  Recommendations 


One  recommendation  stemming  from  the  factors  discussed  above  is  that 
the  Subsurface  Science  Program  should  more  closely  coordinate  microscopic 
and  macroscopic  aspects  of  research.  Specifically,  this  suggestion  is  directed 
at  coordination  of  the  characterization  of  the  microscopic  properties  of  cells 
with  characterization  of  their  adhesion,  and  with  studies  of  mass  transport  in 
both  laboratory  simulation  and  field  studies.  Modeling  of  transport  cannot 
be  carried  out  realistically  without  adhesion  rates  of  bacteria  as  an  input 
parameter,  and  in  converse,  the  significance  of  changes  in  cell  properties  can 
only  be  assessed  ultimately  by  seeing  their  effect  on  transport.  Thus,  we  rec¬ 
ommend  strongly  that  some  microbiological  expertise  be  added  in  the  EMSL 
and  that  some  of  the  resources  of  the  EMSL  be  directed  toward  characteriza¬ 
tion  of  cell  wall  properties  and  their  relationship  to  adhesion.  In  Subsection 
3.3,  some  specific  types  of  characterization  experiments  were  suggested  that 
are  well  within  the  existing  capabihties  of  EMSL.  We  do  not  recommend 
instituting  an  open-ended  fundamental  study  of  bacterial  adhesion  within 
EMSL.  Instead,  we  recommend  a  strongly  focused  collaborative  effort  with 
other  aspects  of  the  subsurface  program,  in  which  adhesion  measurements 
are  made  to  correlate  directly  with  the  bacteria  and  conditions  to  be  used  in 
macroscopic  transport  studies. 
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4  TRANSPORT  OF  BACTERIA  THROUGH 
POROUS  MEDIA 


Many  of  the  proposed  applications  for  bioremediation  involve  the  trans¬ 
port  of  bacteria  through  porous  rock  to  contaminated  zones  underground. 
Transport  through  porous  media  is  an  area  of  quite  active  research  (see  e.g., 
[23])  and  there  are  some  relatively  simple  but  nonetheless  qualitatively  useful 
methods  for  predicting  how  material  will  propagate  under  the  joint  action 
of  advection  by  the  flow  and  diffusion.  In  this  section  we  will  consider  a 
network  model  of  a  porous  rock  formation  and  modify  the  usual  treatment 
of  dispersion  in  order  to  account  for  bacterial  adhesion.  The  results  of  this 
exercise  lead  to  several  serious  concerns  that  must  be  dealt  with  in  any  at¬ 
tempt  at  subsurface  bioremediation;  these  will  be  detailed  at  the  end  of  this 
section. 

Part  of  the  difficulty  in  studying  the  issue  of  bacterial  transport  through 
porous  media  is  that  there  is  no  one  typical  situation.  Instead,  each  subsur¬ 
face  geology,  hydrological  setting,  bacteria  type,  etc.,  leads  to  a  problem 
with  different  considerations  [24].  The  approach  taken  here  is  to  focus  on 
one  possible  situation — transport  of  1  /im  gram  negative  bacteria  through  a 
sandstone  formation  that  is  completely  water  saturated  (i.e.,  below  the  water 
table).  We  will  neglect  bacterial  chemotaxis  (until  the  discussion  section  at 
the  end)  and  assume  a  diffusivity  of  2  x  10“®  cm^/sec.  The  flow  rate  will  be 
taken  to  be  10“®  cm/sec  ~  3  m/year.  If  one  assumes  a  typical  path  length  in 
the  porous  rock  to  be  one  grain  size  (chosen  to  be  2  x  10“^  cm),  the  dimen¬ 
sionless  number  governing  the  relative  importance  of  advection  vs.  diffusion 
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takes  the  form 


Pe  =  —  ~  100. 

This  means  that  advection  is  the  major  means  of  transport  and  dispersion 
(that  is  longitudinal  spreading  of  an  initial  pulse  of  tracer)  is  dominated  by 
the  diffusive  hopping  among  the  various  streamlines  of  the  flow-field.  Figure 
4-1  shows  what  the  concentration  field  of  a  tracer  looks  hke  in  this  high  Peclet 
number  regime. 

To  get  started  towards  the  porous  media  calculation,  we  begin  with 
the  simplest  problem  to  consider,  a  molecular  scale  tracer  in  a  fluid  flowing 
through  a  cylindrical  tube.  The  basic  equation  to  be  solved  is 

{v-V)c  =  D,nV^c,  (4-1) 

where  the  velocity  is  given  by  the  Poiseulle  flow  field 

v  =  2zvm(l-^^.  (4-2) 

Following  Taylor  and  Aris  [25],  one  can  easily  determine  the  effective  diffusion 
constant  as  follows.  Let  us  write 

c  =  c{z)  Sc{r,  z).  (4  —  3) 

We  will  assume  that  the  variation  in  z  (and  <)  is  sufficiently  smooth  to  allow 
the  transverse  profile  6c  to  relax  to  the  solution  of  a  purely  radial  equa¬ 
tion.  Specifically,  we  have  upon  averaging  the  basic  Equation  (4-1)  in  the  r 
direction 

<>- 1  r 


dc 

<  u  >  •Vc(z)-|-  <  V  •  V6c{r,  t)  >  = 
at 

D  — 
^dz^' 

(4-4) 

To  find  Sc  we  use  the  residual  equation. 

(u—  <{?>)•  Vc(5:)  =  Dm  1 

fdHc 

1 

(4-5) 

[  dr^ 

r  dr  j  ^ 
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Scaling  Structure  of  Tracer  Dispersion  Fronts  In... 


Figure  4-1.  (a)  Digitized  photograph  of  the  dispersion  front.  The  experiment  was  conducted  at 
P^clet  number  Pe  ~  4  x  10^.  The  size  of  the  image  is  1000  x  700  pixels. 

(b)  Close-up  photograph  of  the  dispersion  front  showing  concentration  variation 
from  the  sub  pore  size  scale  up  to  the  tens  of  pore  sizes.  White  dots  in  the  photo 
are  glass  beads,  (c)  -  (e)  Development  of  an  equiconcentration  dispersion  front  at 
relative  concentration  C/Co=0.4.  The  fronts  were  taken  at  times:  135  s  (c),  254  s 
(d),  and  400  s  (e).  The  front  at  time  t=400  s  was  reduced  to  a  single  value  function 
by  replacing  parts  of  the  curve  within  overhangs  by  straight  line  segments  in 
jc  direction  (d). 
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where  we  have  dropped  the  z  and  t  derivatives  of  the  residual  profile  in  favor 
of  the  (larger)  radial  derivatives.  We  easily  obtain 

Vm  R^\  dc  , . 


c  _  f’m  (r'^  R^\  dc 


(4-6) 


where  the  last  term  enforces  <  Sc  >=  0.  Note  that  is  automatically 

zero,  i.e.,  there  is  no  net  loss  of  particles  out  of  the  tube.  We  thus  arrive  at 

2D„az^  \  V  3  )/ 

which  then  implies  that  Equation(4-4)  can  be  written  as 

(4-7) 


Deff  —  Dm  + 


vIr:^ 


What  changes  when  we  consider  transport  of  a  bacterium  instead  of  a 
molecule?  There  are  several  different  effects.  One  is  purely  due  to  the  size 
of  the  bacterium.  If  we  take  1/xm  as  a  typical  size,  this  is  hkely  to  be  a  fair 
fraction  of  a  typical  narrow  passage  in  a  porous  rock.  This  complicates  the 
flow  field  and  at  very  least  prevents  the  particle  from  approaching  too  closely 
to  the  wall.  If  we  just  exclude  this  near-wall  region  for  particle  occupation, 
we  derive  to  leading  order  in  \  =  a/ R  (a  =  particle  size) 

/o^"“  v{r)rdr 


/o«-“rdr 


~  U,n(l-H2A). 


Physically,  this  increase  is  due  to  the  lack  of  bacteria  at  the  slowest  moving 
part  of  the  flow  field.  On  the  other  hand,  the  narrowing  of  the  range  of 
streamlines  available  to  the  particle  will  reduce  the  Taylor-Aris  dispersion 
coefficient.  For  small  A,  the  results  have  been  worked  out  by  Brenner  [26] 


with  the  results 


V  —  Um(l  4"  2A  —  4.92A^  ••) 

R2„2 

Deff  =  Dm  +  J^(l-1MX  +  9.68X‘^) 
^oUfn 
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There  are  phenomenological  relations  that  extend  these  small  A  results  to 
A  ~  1;  these  are  quoted  in  Sahimi  [27].  At  least  for  the  medium  we  consider, 
the  average  passageway  will  have  a  radius  of  10  fim,  large  enough  so  that 
these  corrections  are  25%  quantitative  changes.  For  the  smallest  channels, 
the  bacterial  transport  will  be  highly  altered.  However,  these  channels  also 
do  not  support  much  flow  (they  are  like  the  highest  resistances  in  a  random 
resistor  network,  and  the  current  avoids  them)  and  the  advective  dominance 
of  the  transport  means  that  nothing  gets  through  these  pathways  anyway. 
All  told,  we  will  neglect  the  finite  size  constraint  for  the  case  at  hand.  This 
will  most  certainly  not  be  justified  in  situations  with  less  permeable  rocks 
and/or  much  slower  flow  rates. 

The  next  issue  to  consider  is  the  existence  of  nontrivial  electrostatic 
forces  associated  with  (negative)  charges  on  the  surface  of  the  cell.  These 
forces  are  strongly  dependent  on  the  ionic  strength  of  the  aqueous  environ¬ 
ment  and  also  the  nature  of  the  various  surfaces  involved;  these  have  been 
discussed  at  length  in  other  sections  of  this  report.  One  way  of  seeing  this 
effect  is  to  consider  the  change  in  average  velocity  due  to  finite  size  effects, 
but  now  including  a  charge  density  of  <t  ~  10^^  charges/cm^  [28].  At  solu¬ 
tions  with  less  than  10“^  M  electrolyte  concentration,  we  see  in  Figure  4-2 
the  expected  behavior  of  a  slight  increase  in  average  velocity  as  the  size  of 
the  bacterium  is  increased,  due  to  the  aforementioned  exclusion  from  the 
slowest  streamlines.  As  the  molarity  is  increased  to  numbers  typical  of  sea 
water,  there  is  a  change  to  the  opposite  behavior.  This  is  simply  due  to  the 
trapping  of  particles  at  the  secondary  Van  der  Waals  minimum. 

Typical  groundwater  has  an  electrolyte  concentration  in  the  range  10“^  - 
10“^  M.  We  will  assume  that  there  is  no  strong  minimum  and  hence  no  purely 
electrostatic  (reversible)  adsorption.  Our  reason  for  this  is  quite  simple — if 
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Figure  4-2.  Dependence  of  the  velocity  ratio  on  particle  radius  and  the  electrolyte  concentration 
of  the  carrier  solution.  Ratios  less  than  unity  result  from  nonuniform  sampling  of 
radial  positions  by  the  Brownian  particle  in  favor  of  the  slow-moving  region  near  the 
capillary  wall.  Nonuniform  sampling  is  caused  by  the  lower  potential  energy  of 
interaction  between  particle  and  capillary  at  separations  corresponding  to  the 
secondary  minimum. 
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there  is  such  an  adsorption,  all  the  bacteria  would  reach  the  wall  in  about  the 
diffusion  time  for  a  single  pore  and  everything  would  stick  irreversibly  on  a 
time  scale  very  short  compared  to  the  needed  transit  time  (hours  vs.  months). 
(Indeed  this  is  what  is  revealed  by  Table  4.1  which  shows  the  transition  from 
reversible  to  irreversible  adsorption  of  a  Pseudomonas  strain  occurring  in 
about  one  day.)  In  other  words,  the  only  hope  for  transport  is  for  adhesion 
to  be  a  very  rare  event  which  is  not  mediated  by  strong  physadsorption.  We 
will  therefore  proceed  with  the  assumption  that  we  can  neglect  electrostatic 
forces. 

Another  possible  complication  is  that  pore  walls  are  quite  rough  and 
can  create  small  eddies  in  the  flow  as  well  as  snag  bacteria  on  the  high 
points  [29].  It  is  quite  difficult  to  take  into  account  this  roughness  as  part 
of  the  hydrodynamic  calculation.  On  the  other  hand,  the  actual  nature  of 
the  surface  will  automatically  be  taken  care  of  in  the  sticking  coefficient  (see 
below)  if  measurements  are  made  with  the  relevant  geological  materials. 

So,  we  are  left  with  the  one  predominant  effect  that  one  must  consider 
when  trying  to  estimate  the  feasibility  of  bioremediation.  As  the  bacteria 
transit  through  the  porous  media,  some  of  them  will  irreversibly  adhere  to 
the  pore  walls.  This  may  happen  whether  or  not  the  bacteria  are  in  a  dormant 
phase — in  fact,  there  are  cases  when  the  dormant  phase  is  a  better  adherent. 
In  regions  that  lie  between  the  input  and  the  contaminated  zone,  adsorption 
represents  a  net  loss  of  useful  biomass.  In  contaminated  regions,  adsorption 
affects  the  spatial  distribution  of  bacterial  deposition  and  the  range  over 
which  we  can  use  advection  to  scatter  the  bacteria  and  not  rely  on  the  much 
slower  process  of  contaminant  diffusion.  Let  us  assume  for  simplicity  that  we 
are  never  in  the  situation  of  having  a  large  fraction  of  the  possible  sticking 
sites  filled  with  bacteria.  Then  the  local  flux  through  the  surface  can  be 
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Tabic  4-1. 


Proportions  of  Pseudomonas  Strain  EK20  Cells  Adhering  Reversibly  and  Irreversibly^ 
to  a  Polystyrene  Substratum  at  Different  Time  Intervals 


%  of  total  bacteria  attached 


Time  (h)  Reversibly  Irreversibly 


89.5 

10.5 

70.7 

29.5 

44.2 

55.8 

27.3 

72.7 

22.0 

78.0 

24.0 

75.6 

9.4 

90.6 

a  Determined  by  using  3H-labeled  bacterial  suspensions  in  polystyrene  petri  dishes.  After  the 
supernatant  was  poured  off,  the  washings  were  retained  as  the  reversibly  adhered  fraction  while 
the  irreversibly  adhered  fraction  remained  attached  to  the  substratum. 


taken  to  be 


=  adoc, 
or 


where  do  is  the  thickness  of  the  layer  over  which  adhesion  could  occur  (pos¬ 
sibly  a  few  tens  of  nm)  and  a  is  a  rate  constant. 

We  now  redo  the  above  Taylor-Aris  calculation.  The  equation  for  6c 
now  has  the  more  general  solution 


Using  the  above  boundary  condition  yields 


dc  _  D  fB\ 
8D  dz  ado\R)' 


The  other  condition  that  needs  to  be  enforced  is 


<  6c  >=  0, 


which  yields 


Solving,  we  obtain 


2^  12D  dz 


n  —  1  VmR^  dc 
^  12  D  dz 

^  ^  adoR 
VmR^  dc  / 


24D  dz 
2D 


/o  \ 

V  adoIU  ■ 


The  critical  parameter  which  governs  this  differs  from  the  previous  case  is 
2DfadoR.  If  we  substitute  the  typical  numbers  we  have  been  using,  we  get 
that  this  parameter  is  numerically  equal  approximately  to  (if  we  measure 
time  in  units  of  seconds).  Since  a~^  is  fairly  large,  we  can  approximate  the 
part  of  B  that  depends  on  c  as 

—adoR 
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Now  the  B  term  in  8c  represents  a  net  outflow  of  material  from  the  mean  flow. 
That  is,  the  average  value  of  ^  longer  zero  because  of 

the  8  function  arising  from  the  In  r  term  in  8c.  It  is  easy  to  see  that  the  net 
flow  equation  becomes 

dc  .  dc  d'^c 
dz^ 

where  there  is  a  slight  modification  to  and  a  net  loss  term 


Note  that  the  diffusion  constant  has  dropped  out,  since  the  loss  is  reaction- 
limited.  In  the  limit  of  large  a  (for  example  if  there  is  strong  physadsorption 
due  to  large  electrolyte  content  in  the  water),  we  would  have  obtained  the 
diffusion-limited  result  ^  ~  SirDmlR^-  For  our  case,  we  are  using  do  ~  10“® 
cm,  R  ~  10“^  cm,  and  a  ~  (lO"**  —  10"^),  giving  ^  in  the  range  10“®. 


We  now  put  this  single-channel  result  into  a  network  model  of  a  porous 
medium.  We  use  a  specific  model  due  to  Koplik  et  al.  [30]  where  the  flow  is 
along  bonds  on  a  square  (in  three  dimensions  a  cubic)  lattice.  These  bonds 
are  “allowed”  with  probability  p  and  there  is  a  critical  value  of  p  (the  bond 
percolation  threshold)  above  which  there  is  a  cluster  that  continuously  con¬ 
nects  an  input  node  to  an  output  node.  The  widths  of  the  bonds  are  chosen 
from  a  uniform  distribution;  we  have  chosen  a  mean  channel  size  of  lO/^m 
with  the  distribution  running  to  5pm  to  15/im.  The  calculation  proceeds  as 
follows.  First,  the  permeability  of  each  channel  is  determined  from  its  width 
and  the  flow  problem  is  solved  via  an  analogous  resistor  network  (as  already 
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mentioned,  pressure  drop  is  like  voltage,  fluid  flux  like  current).  This  deter¬ 
mines  the  velocity  in  each  bond  as  well  as  the  ratio  of  fluid  flux  to  pressure 
drop  for  the  entire  porous  medium  model.  In  physical  units,  the  perme¬ 
ability  of  the  medium  equals  ^^/(p)  where  R  is  the  average  channel  radius 
{R  =  10pm),  d  the  channel  length  (200pm  here)  and  f{p)  is  a  dimensionless 
function  calculated  by  the  program.  Substituting,  we  get 

A;  =  9.4  X  10”^°  cm^  f{p) 

and  f{p)  for  a  20  x  20  lattice  is  shown  in  Figure  4-3.  For  comparison, 
permeability  in  sandstone  ranges  from  10“®  to  10"^^.  For  definitions,  we  will 
use  p  ~  .7  in  our  simulation. 


Given  the  velocities  in  each  channel,  we  can  now  solve  the  mean  con¬ 
centration  transport  equation  to  determine  how  the  bacteria  move  through 
the  medium.  The  easiest  way  to  do  that  is  to  take  the  Laplace  transform  of 
the  concentration  equations  to  reduce  everything  to  a  matrix  problem  which 
depends  on  the  transform  variable  s.  To  see  how  this  works,  let  us  consider 
a  channel  that  connects  nodes  i  and  j.  The  equation  to  be  solved  is 


dc  ^  ^ 


The  solution  takes  the  form 


c(z;  s) 
a± 


u  ±  -I-  4D{s  + 

W 


The  boundary  conditions  involve  fixing  the  nodal  concentration  c(0,  s)  = 
Cj;  c(/,  s)  =  Cj.  After  some  algebra,  we  find  that  one  can  express  the  flux 
from  i  to  j  as 
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Figure  4-3.  Dimensionless  permabiliiy  from  bond  percolation  network  model. 
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(where  Sij  is  the  cross  sectional  area)  and  where 

DSi 

Gfj  =  — (m  +  nig  coth  m^) 
nc.. 

Gfj  =  — (mje"”*  sinh  m,) 


m  =  vll2D  (the  local  Peclet  number) 
m,  =  (m"  +  (s  +  y8)/Vl>)i/^ 

Given  a  source  at  some  overall  input  node  i  =  0  one  can  solve  the  set 
of  network  equations 

=  E  Jir 

} 

The  left-hand  side  is  the  Laplace  transform  of  a  ^  function  (in  time)  at  the 
input  and  all  non-input  nodes  cannot  accumulate  nor  create  material.  At  the 
output  node  N,  the  flux  —  J^j  represents  the  first  passage  time  distribu¬ 
tion.  One  can  either  use  an  inverse  Laplace  transform  method  to  recover  P{i) 
or  more  simply,  calculate  moments  of  P{t)  via  the  small  s  expansion  of  its 
Laplace  transform.  The  only  change  that  we  had  to  implement  (as  compared 
to  the  original  use  of  this  model  for  tracer  dispersion)  was  the  introduction 
of  non-zero  decay  rate  p.  It  is  easy  to  show  that  for  =  0,  f  P(t)dt  =  1 
corresponding  to  /*(s)|s=o  =  1-  Hence,  the  net  loss  of  material  as  it  passes 
through  the  porous  medium  is  precisely  1  —  P{s  =  0).  Similarly  one  can 
obtain  expressions  for  the  mean  passage  time 


<  t  >  -- 


j  tp(t)i  I  p{t) 


as  well  as  a  measure  of  the  longitudinal  dispersivity  <  P  >  —  <  t 


For  the  case  we  have  used  as  our  model  calculations,  we  obtained  a 
throughput  of  .89  in  a  20  x  20  lattice.  Since  the  bond  length  is  200/xm, 
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this  corresponds  to  roughly  a  10%  loss  in  4  m  of  rock.  The  average  time  of 
passage  through  the  material  was  4  x  10^  seconds,  with  a  standard  deviation 
due  to  dispersion  of  1  x  10^  seconds.  These  numbers  are  fairly  alarming.  If 
one  wants  to  go  through  4  mm  of  rock,  one  quickly  estimates  that  only  1  in 
e+ioo  bacteria  survive  the  journey. 

What  implications  does  this  have  for  subsurface  remediation?  The  first 
point  is  that  estimates  based  on  simple  conceptualizations  of  porous  media 
are  probably  more  useful  at  this  stage  than  full-blown  supercomputer  cal¬ 
culations.  Models  of  this  sort  can  be  used  to  study  the  distribution  (and 
the  sticking  locations)  of  bacteria,  co-transported  nutrients  and  in  fact  the 
original  contaminants.  For  example,  the  fact  that  bacteria  are  transported 
at  large  Peclet  number  but  small  molecules  with  diffusivities  of  10“®  cm^/sec 
are  transported  at  small  Peclet  numbers  immediately  leads  to  a  different  dis¬ 
tribution  inside  the  rock  (the  small  molecules  will  explore  much  more  of  the 
“dead-end”  parts  of  the  pore  space  that  are  not  on  the  main  flow  path — they 
will  therefore  suffer  much  more  dispersion  as  they  move  through  the  porous 
material. 

Specifically  we  feel  that  the  most  glaring  uncertainties  lie  in  the  stack¬ 
ing  coefficient  and  in  fact  in  the  entire  realm  of  how  much  of  the  microscopic 
behavior  is  determined  by  relatively  unchanging  physics  and  chemistry  vs. 
being  determined  by  microbiology,  which  can  make  vast  changes  in  cell  phys¬ 
iology,  (under  the  time  scales  we  are  considering  [31]).  Therefore,  one  should 
at  this  stage  be  using  simple  porous  medium  models  and  estimates  (and  sim¬ 
ilarly  simple  experimental  realizations  of  porous  media  such  as  micromodels 
[32]  or  bead  packs)  instead  of  trying  to  do  large-scale  computations  and/or 
experiments  in  complex  geomorphologies,  and  having  the  results  totally  de¬ 
pendent  on  unknown  microbiological  issues. 
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The  second  point  is  that  we  feel  the  issue  of  irreversible  adhesion  is 
sufficiently  serious  that  it  might  preclude  the  use  of  bioremediation  in  many 
cases.  One  should  therefore  develop  ideas  for  limiting  or  at  best  eliminat¬ 
ing  this  effect  [33].  One  direct  possibility  is  drilling  lots  of  holes,  ensuring 
direct  access  (within  several  meters  at  worst)  to  contaminated  sites.  This 
is  clearly  expensive.  One  could  help  alleviate  loss  if  one  drills  several  wells 
and  used  high  pressure  to  generate  much  higher  than  ambient  flow  velocities 
that  would  distribute  the  bacteria  more  widely,  all  else  remaining  the  same. 
Pumping  into  rock  formations  is  a  standard  oil  industry  technology  used  dur¬ 
ing  secondary  oil  recovery.  Since  the  flow  rate  would  only  be  enhanced  to  a 
large  degree  near  the  well,  many  holes  would  again  be  needed. 

Another  possible  mitigation  approach  would  be  encapsulation  of  the 
bacteria  in  non-sticky  containers.  In  some  sense,  this  is  a  chemical  means  of 
getting  a  non-adhesive  coating  on  a  particular  bacterium.  One  might  also 
imagine  doing  the  same  thing  by  genetically  engineering  a  bacterium  and 
thereby  changing  its  sheath  of  glycolipids  (and  other  stuff).  For  the  capsule 
idea,  there  is  a  stringent  size  constraint.  We  are  already  pushing  the  limits 
of  size  in  getting  bacteria  through  the  pore  space  of  sedimentary  rock  (frac¬ 
tures  are  probably  wider  but  at  least  in  horizontal  transport  these  are  not 
the  pathways  that  dominate  contaminant  transport — hence  the  bacteria  go 
to  the  wrong  spots).  Adding  a  coating  (say,  placing  a  bacterium  in  a  2  /xm 
radius  hollow  hposome — ^a  self- assembled  object  made  from  pure  phospho¬ 
lipids)  makes  matters  worse.  Using  dormant  phases  of  the  bacteria  would 
help  (these  are  usually  significantly  smaller)  but  then  one  must  ensure  that 
they  are  “revived”  upon  reaching  the  right  spot. 

The  next  point  is  a  statement  that  follows  from  the  dissimilarity  in  con¬ 
taminant  distribution  as  compared  to  any  conceivable  bacterial  distribution 
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obtained  by  physical  means.  Because  the  contaminants  have  explored  tens  of 
years  worth  of  possible  paths  in  the  porous  rocks  and  because  they  are  much 
smaller  than  the  bacteria  (as  already  mentioned  this  greatly  affects  their 
transport),  even  under  the  best  of  circumstances  one  will  have  to  wait  for 
contaminants  to  diffuse  out  of  narrow  dead-ends  (and  possibly  to  de-adhere) 
under  the  action  of  concentration  gradients  created  by  having  the  bacteria 
degrade  all  the  contaminant  in  the  bacteria- accessible  regions.  This  means 
that  it  takes  at  least  as  long  to  clean  up  a  site  as  it  did  to  pollute  it  in  the  first 
place.  There  can  be  no  magic  potion  based  on  diffusion  of  the  remediator 
that  undoes  years  of  diffusion  of  the  contaminant  without  etching  the  rock 
to  change  the  pore  structure. 

The  final  point  is  probably  the  most  general.  That  is,  we  have  in  this 
chapter  been  treating  bacteria  as  “living  colloids”  with  some  fixed  set  of 
physical  and  chemical  properties.  On  a  short  time  scale  (say  hours)  this  is 
quite  reasonable.  As  we  start  talking  about  processes  that  will  last  years 
to  tens  of  years,  microbiology  becomes  more  relevant  than  physics.  This  is 
both  a  blessing  and  a  curse.  One  might  be  able  to  rely  on  colony  growth 
by  division,  active  de-adhesion  (if  the  food  supply  goes  down,  bacteria  will 
actively  disassociate  themselves  from  walls  to  try  to  find  more  congenial 
surroundings),  chemotaxis,  cell  adaption  to  local  environments,  etc.  If  the 
cell  “wants”  to  find  the  contaminant  (i.e.,  it  makes  beneficial  use  of  the 
contaminant  in  some  sort  of  metaboHc  process),  it  may  succeed  even  if  the 
forces  of  physics  and  chemistry  are  somewhat  predisposed  against  it.  It  is  of 
course  a  major  challenge  to  figure  out  how  to  get  bacteria  “interested”  in  a 
particular  undesirable  (from  the  human  perspective)  material.  And  this  leads 
directly  to  the  curse  part  of  the  dominance  of  microbiology.  On  the  same 
timescale,  any  bioremediant  will  have  to  compete  in  the  particular  subsurface 
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ecology  with  all  of  the  native  organisms.  Just  designing  (or  finding)  the 
right  organism  will  not  do  —  something  which  happily  digests  EDTA  but 
cannot  compete  with  the  natives  in  the  particular  aquifer  which  has  EDTA- 
heavy  metal  contamination  is  not  useful  for  bioremediation  purposes.  This 
is  a  particularly  difficult  challenge  since  our  understanding  of  the  microbial 
ecology  of  the  subsurface  is  essentially  nil.  The  only  logical  way  to  proceed  is 
by  working  with  bacteria  already  adapted  for  the  particular  environment  and 
“breed”  into  them  (either  by  genetic  engineering  or  by  relying  on  selection  in 
the  proper  contaminated  environment)  the  ability  to  degrade  the  offending 
“stuff” .  One  then  has  to  hope  that  utilizing  this  stuff  gives  them  enough  of 
a  competitive  advantage  actually  to  live  long,  prosper,  and  thereby  clean. 
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5  METAL  BIOPROCESSING 


Depending  on  context,  metals  are  hazardous  wastes  to  be  accumulated 
and  disposed  of,  or  economic  assets  to  be  recovered  from  waste  or  extracted 
from  ores.  Microorganisms  have  a  role  to  play  in  both  of  these  tasks  and  have 
attracted  considerable  attention  in  this  connection  over  the  recent  decades. 
Metal  bioprocessing  has  been  applied  to  the  remediation  of  contaminated 
surface  and  ground  waters,  to  the  treatment  of  industrial  waste  streams,  and 
in  the  extractive  and  waste  treatment  operations  of  the  mining  industry  [34]. 

For  bioremediation  of  waters  containing  solubilized  heavy  metals,  the 
various  processes  of  interest  in  some  cases  involve  metabolically  active  mi¬ 
crobes,  in  other  cases  non-living  (i.e.,  non-metabolic)  organisms.  In  precip¬ 
itation  processes,  soluble  metals  interact  with  substances  secreted  by  living 
microbes  to  produce  insoluble  metal  compounds.  Sulfate-reducing  bacteria 
such  as  Desulfovibrio  and  Desulfotomaculum  are  a  prominent  example.  In 
anoxic  environments  with  an  abundance  of  sulfate  and  degradable  organic 
matter,  they  induce  the  production  of  H2S  and  bicarbonate: 

2CH2O  -I-  SO^  ^  H2S  -h  2HC0_3.  (5  -  1) 

The  sulfide  interacts  with  soluble  divalent  metal  ions  (M)  to  produce  insoluble 
metal  sulfides: 

M^"*"  -f  MS  (precipitate).  (5  —  2) 

These  processes  are  relevant  to  the  treatment  of  metal-contaminated  waste 
waters  with  high  concentrations  of  sulfate  and  organic  matter,  and  in  some 
places  occur  naturally.  In  many  applications,  e.g.,  to  the  acidic  drainages 
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from  mining  operations,  it  may  be  necessary  to  add  the  necessary  organic 
compost  deliberately  or  in  other  ways  induce  the  growth  of  organic  matter. 

Intracellular  accumulation  refers  to  processes  in  which  soluble  metals 
first  attach  to  the  cellular  surface  of  the  microorganism  and  are  then  trans¬ 
ported  across  the  cell  membrane  to  the  interior,  where  they  may  precipitate 
or  bind  to  cellular  components.  Transport  across  the  cell  membrane  in  most 
cases  is  thought  to  involve  energy-dependent  processes  in  hving  organisms. 
These  processes  are  not  yet  well  understood.  It  may  be  that  the  heavy  met¬ 
als  are  able  to  confuse  the  cell  by  using  the  same  mechanisms  that  transport 
metaboUcally  essential  ions  such  as  potassium.  The  uptake  of  heavy  met¬ 
als  by  this  process  is  sensitive  to  a  number  of  environmental  factors  in  the 
wastewater,  e.g.,  pH,  the  presence  of  molecules  that  complex  the  metals,  etc. 
Substantial  intracellular  accumulations  have  been  observed.  For  example, 
with  a  mixed  culture  consisting  of  Pseudomonas  maltophilia,  Staphylococcus 
aureus  and  a  cornyeform  organism,  accumulation  of  Ag"*"  to  a  level  of  30%  of 
biomass  dry  weight  has  been  observed. 

Metal  immobilization  through  oxidation  and  reduction  reactions  can  be 
illustrated  by  the  example  of  oxidation  of  manganese  at  chemically  reactive 
sites  on  the  cell  walls  of  Bacillus  spores: 

Mn=*+  -I- 1/2  O2  +  H2O  Mn02  (insoluble)  -|-  2H+.  (5  -  3) 

Exemplifying  a  reduction  process,  the  reduction  of  gold  from  Au^"^  to  Au® 
has  been  detected  on  the  cell  walls  of  Bacilli  subtilis.  Extracellular  com- 
plexation  of  metals  is  another  mechanism  by  which  bacteria  immobilize  sol¬ 
uble  metals.  This  arises  where  the  organisms  excrete  appropriate  chelating 
agents,  for  example  the  siderophores  that  can  be  isolated  from  Pseudomonas, 
Arthrobacter  and  other  bacterial  species.  In  some  cases,  there  is  the  prospect 
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that  the  complexing  agents  can  be  recovered  and  reused  repeatedly  for  metal 
accumulation. 

The  attachment  of  soluble  metals  to  microbial  cells  or  to  the  substances 
excreted  by  cells  is  typically  quite  sensitive  to  the  composition  and  chemical 
state  of  the  wastewater  that  is  being  treated.  Organic  and  inorganic  con¬ 
taminants  can  affect  the  chemistry  and  availability  of  the  metals  for  binding 
to  the  bioabsorbents.  The  latter  may  be  subject  to  attack  by  chemical  sub¬ 
stances  in  the  solution,  etc.  This  implies  that  the  appropriate  technology  for 
bioremediation  of  metals  may  have  to  be  finely  tailored  to  the  particulars  of 
the  varying  applications. 

A  number  of  bioremediation  processes  have  been  developed  and  intro¬ 
duced  for  commercial  applications.  One  of  these  is  the  BIOCLAIM  process. 
This  is  based  on  microorganisms,  principally  Bacillus,  whose  cells  are  de¬ 
liberately  ruptured  in  a  strong  caustic  solution  to  expose  functional  groups 
with  good  metal  binding  capacities.  The  cell  fragments  are  immobilized  on 
extruded  polyethyleneimine  beads.  Other  processes,  similarly  involving  non¬ 
living  mixtures  of  biomass  that  are  immobilized  on  beads  that  can  be  used 
in  fixed-bed,  fluidized-bed  or  dispersed-bed  reactors,  go  under  the  names  Al- 
gaSORB  and  BIO- FIX.  At  the  present  time,  most  of  the  systems  that  are 
in  commercial  use  or  that  are  approaching  commercialization  are  based  on 
non-living  bacteria.  Living  systems  have  the  advantage  that  they  are  a  re¬ 
newable  resource  and  that,  moreover,  they  exploit  metabolic  processes  (e.g., 
the  reduction  of  sulfate  to  sulfide)  not  mediated  by  non-living  systems.  The 
difficulty  is  that  waste  streams  are  often  toxic  to  hving  systems  and  devoid 
of  nutrients. 

Biotechnology  is  likely  to  have  an  increasingly  important  role  in  mining. 
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an  industry  in  the  U.S.  that  is  in  the  course  of  long-needed  modernization. 
Exploitation  of  increasingly  lower-grade  and  refractory  ores  is  one  of  the 
big  challenges  to  the  industry.  The  refractory  ore  problem  is  well  exem¬ 
plified  by  the  case  of  gold  mining.  Oxide  gold  ores,  which  have  been  the 
backbone  of  the  industry,  respond  well  to  the  standard  cyanide  method  of 
extraction.  Here  the  gold  component  in  crushed  ores  complexes  with  cyanide 
and  the  complex  is  removed  by  absorption  on  activated  carbon.  But  reserves 
of  oxide  ores  are  being  depleted  worldwide  and  are  being  replaced  by  ores 
that  consist  of  sub-micron  particles  of  elementary  gold  embedded  in  a  sul¬ 
fide  mineral  matrix,  e.g.,  pyrite  (FeS2).  The  cyanide  does  not  penetrate  well 
here  and  extractive  yields  are  small,  often  less  than  20%  of  the  available 
gold.  Alternative,  non-biological,  technologies  exist,  but  they  are  costly.  A 
biological  alternative — bioleaching — involves  the  use  of  aerobic,  acidophilic 
bacteria,  e.g.,  Thiobacillus,  which  oxidize  sulfur  compounds,  degrading  the 
sulfide  matrix  and  exposing  the  elemental  gold.  Bioleaching  of  refractory 
gold  concentrates  is  being  pursued  commercially  in  South  Africa,  Australia, 
and  South  America.  At  the  present  time,  the  largest  application  of  biotech¬ 
nology  in  mining  is  to  the  bioleaching  of  low-grade  sulfidic  copper  ores.  This 
is  now  practiced  widely  throughout  the  western  United  States  and  especially 
in  Chile,  chiefly  for  ores  containing  chalcocite  (CU2S)  and/or  covellite  (CuS). 

Tightened  restrictions  on  sulfur  dioxide  emissions  from  coal-fired  power 
plants  have  increased  attention  to  pre-combustion  desulfurization  technol¬ 
ogy.  There  would  appear  to  be  a  vast  market  for  coal  biodesulfurization  if 
competitive,  cost-effective  techniques  could  be  developed.  Research  with  the 
thermophilic,  acidophilic  bacterium  Sulflobus  acidocaldarius  for  extraction 
of  sulfide  and  organically  bound  sulfur  has  been  performed.  At  the  present 
time,  however,  commercially  promising  findings  have  not  emerged. 
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6  IN  SITU  RECOMBINANT  DNA 


In  order  to  accomplish  in  situ  bioremediation  one  often  is  faced  with 
the  need  specifically  to  “engineer”  an  organism  that  will  have  the  desired 
metabolic  behavior  for  the  contaminants  at  hand.  Normally,  this  engineering 
is  imagined  to  occur  in  the  laboratory,  creating  an  organism  that  is  both 
adapted  to  living  in  the  specific  in  situ  environment  (in  terms  of  temperature, 
pressure,  chemical  composition,  etc.)  and  able  to  accomplish  the  desired 
degradation.  In  this  section,  we  explore  a  speculative  but  somewhat  novel 
concept — ^that  of  using  the  resident  microorganisms  and  transmitting  only 
the  new  genetic  instructions  from  the  surface. 

The  first  question  to  be  addressed  concerns  the  possible  benefits  such 
a  scheme  might  have  over  more  traditional  concepts.  Essentially,  it  may 
be  quite  difficult  to  transport  desired  microorganisms  from  the  laboratory 
to  a  specific  subsurface  environment.  This  issue  is  discussed  earlier  in  this 
report  in  more  quantitative  detail,  but  the  qualitative  concerns  that  enter 
are:  (a)  will  the  microorganism  be  able  to  survive  transit  through  unfavorable 
regions  to  reach  contaminated  zones;  (b)  will  the  microorganism  adhere  to 
the  surface  of  the  porous  media  to  such  an  extent  as  to  make  long  distance 
transport  unattainable;  (c)  will  heterogeneity  of  a  contaminated  zone  on  the 
scale  of  meters  make  predictions  of  organism  viability,  made  via  tests  on  core 
samples,  rather  unreliable;  and  finally  (d)  will  hard-to-reach  spots  in  the 
porous  medium  remain  locations  of  long-term  contamination  simply  because 
microorganisms  can  no  longer  fit  through  the  small  throats  of  diagenetically 
altered  aquifers. 
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The  first  two  points  of  this  list  might  be  addressed  by  encapsulating 
the  bacteria  for  the  journey  to  the  relevant  location.  This  idea,  coupled 
with  the  concept  that  one  will  quite  often  drill  wells  directly  through  the 
rock  formation  and  pump  the  relevant  bioremedial  “stuff”  directly  into  the 
correct  stratum,  should  certainly  work  in  some  cases.  The  basic  problem, 
though,  is  that  even  a  fairly  tightly  fitting  capsule  of,  say,  5  fim  diameter 
will  have  trouble  negotiating  through  the  pore  space  of  a  typical  aquifer.  So, 
though  one  might  diminish  wall  adhesion  by  proper  capsule  choice,  one  might 
just  get  stuck. 

On  the  other  hand,  sending  encapsulated  genetic  information  solves 
problems  (a)  and  (b)  without  imposing  this  severe  size  constraint.  The  ge¬ 
netic  information  would  be  contained  inside  a  vector,  either  a  virus  known  to 
infect  some  specific  microorganism  or  a  plasmid;  these  could  easily  fit  inside 
capsules  of  size  .1  /xm.  The  best  choice  for  capsules  would  be  self- assembled 
liposomes,  which  could  be  induced  to  form  spontaneously  in  a  “soup”  of  ge¬ 
netic  material  that  would  be  straightforwardly  incorporated  into  the  lipid 
vescicle.  These  can  easily  be  made  in  this  size  scale. 

Finally,  use  of  already  existent  microorganisms  will  help  alleviate  (c) 
(since  one  could  use  the  locally  adapted  colonies)  and  possibly  even  (d)  if 
the  colonies  predate  the  diagenetic  activity  that  led  to  the  “dead  ends”  of  the 
pore  space.  Of  course,  solving  all  the  problems  of  using  laboratory  generated 
organisms  by  shifting  to  “reprogramming”  in  situ  colonies  merely  transfers 
the  difficulty  to  getting  the  right  instructions  to  the  right  place. 

Let  us  assume  that  we  will  be  using  “naked”  DNA  in  the  form  of  a 
plasmid  to  do  the  requisite  genetic  engineering.  We  need  to  find  a  plasmid 
that  is  likely  to  be  able  to  function  in  the  variety  of  related  species  to  be  found 
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in  situ.  These  do  exist — for  example,  the  famous  R  plasmid,  which  confers 
antibiotic  resistance,  will  function  in  a  wide  variety  of  gram-negative  bacteria. 
Next,  we  must  ensure  that  cells  that  incorporate  the  new  genetic  information 
are  “selected”  for  by  the  environment  in  such  a  way  as  to  compensate  for 
the  disadvantage  of  carrying  extra  DNA.  The  obvious  holy  grail  is  that  the 
plasmid  will  allow  the  modified  bacteria  to  grow  more  rapidly  by  specifically 
allowing  them  to  metabolize  the  contaminant  that  one  wants  to  degrade.  If 
this  could  be  accomplished,  one  would  modify  the  ambient  ecology  only  in 
the  contaminated  zones,  and  then  the  modification  itself  would  be  selected 
against  as  the  environment  becomes  decontaminated. 

What  are  the  hard  parts?  Beside  finding  the  relevant  gene  insert  and 
vector  (these  would  have  to  be  done  in  any  genetic  engineering  solution), 
the  actual  incorporation  might  be  tricky.  It  is  known  how  to  do  this  under 
controlled  laboratory  conditions,  and,  in  fact,  the  incorporation  of  DNA  into 
bacteria  was  the  original  experimental  finding  of  Avery  that  identified  DNA 
as  the  carrier  of  genetic  information.  Whether  this  could  be  done  with  large 
enough  yield  in  an  uncontrolled  setting  is  probably  unknown.  One  could  co¬ 
transport  any  specifically  needed  chemicals  to  make  the  bacteria  receptive  to 
naked  DNA,  but  one  would  have  to  make  the  right  mix  in  a  unpredictable 
chemical  bath.  There  is  certainly  no  difficulty  in  testing  this  experimentally 
using  known  host- vector  systems.  One  would  also  have  to  design  the  encap¬ 
sulated  soup  to  allow  for  long-term  (weeks  to  months)  survival  of  the  plasmid 
during  the  transit  phase  of  the  remediation  effort. 

The  difficulties  in  getting  this  idea  to  work  are  non-negligible.  One 
should  keep  in  mind,  though,  that  natural  genetic  engineering  does  occur, 
and  cells  are  quite  happy  to  pick  up  bits  of  useful  information  if  it  helps 
them  to  cope  better  with  their  environment.  What  this  probably  means  is 


71 


that  even  though  indigenous  microorganisms  could  eventually  figure  out  on 
their  own  how  to  “eat”  a  particular  contaminant,  usefully  speeding  up  this 
process  by  directly  providing  the  relevant  information  should  be  attainable 
if  the  actual  competitive  advantage  to  the  organism  is  large  enough  so  that  a 
very  small  yield  will  nonetheless  lead  eventually  to  a  large  number  of  modified 
bacteria. 
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7  THE  STATISTICS  OF  CENSUS  TAKING 


How  is  one  to  estimate  the  subsurface  population  of  the  earth? 

Well,  how  does  the  Census  Bureau  conduct  a  census — (poorly) — but 
apart  from  that  how  do  they  go  about  it?  They  go  to  homes,  and  count 
residents.  (Thus  they  miss  the  homeless,  which  is  hard  on  New  York  State 
when  it  comes  to  political  desserts.)  Even  so,  they  miss  some  residents 
because  no  one  was  at  home  when  they  called — they  miss  others  because 
some  people  do  not  want  to  be  counted  or  identified  in  any  way. 

The  problem  of  taking  a  microbial  subsurface  census  faces  many  of  the 
same  difficulties,  but  we  may  start  by  finding  the  homes.  To  do  this  we 
need  a  reasonably  accurate  description  of  the  subsurface  environment  world¬ 
wide,  and  to  a  fairly  great  depth  (kilometers).  This  description  must  include 
availabihty  of  water,  porosity,  availability  of  assorted  minerals,  etc. 

Once  this  daunting  task  is  completed  (it  includes  the  tundra,  too,  of 
course,  as  well  as  the  subsoil  under  the  Antarctic  ice  sheets),  we  can  roughly 
estimate  the  microbial  population  on  the  basis  of  two  considerations: 


1.  In  an  environment  which  is  roughly  comparable  to  one  we  have  seen 
before  and  for  which  we  know  microbial  density,  the  new  environment 
will  have  microbial  residency  much  like  the  old. 

2.  In  any  environment  where  there  is  a  possibility  of  migration  due  to 
porosity  or  flow,  and  in  which  there  is  available  chemical  energy  in  the 
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mineral  landscape,  one  will  find  substantial  density  of  microbial  life 
exploiting  the  available  energy. 

Occasional  sampling  of  the  subsurface,  of  necessity  a  very  thin  sampling, 
will  serve  to  reinforce  or  deny  one’s  estimates  based  on  principle  ^2.  Sam¬ 
pling  itself  will  never,  we  should  think,  give  enough  data  to  estimate  reliably 
the  population  statistics,  in  the  absence  of  a  very  strong  initial  hypothesis. 
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8  NUCLEAR  WASTE  DISPOSAL— BOROSILICATE 
GLASS 

8.1  Introduction 


U.S.  nuclear  waste  currently  in  temporary  repositories  will  be  incorpo¬ 
rated  into  borosilicate  glass  (BSG)  for  ultimate  storage  in  a  geologic  repos¬ 
itory.  BSG  was  chosen  for  solidification  of  high-level  radioactive  waste  at 
the  Savannah  River  Plant  (SRP)  over  a  number  of  alternative  glass  (e.g., 
phosphate  glass),  glass  ceramic,  ceramic  (e.g.,  synrock,  other  ceramics  tai¬ 
lored  specifically  for  the  waste  stream),  and  mixed  materials  developed  in 
the  1970s  and  1980s  [35].  The  Australian  waste  form  “synrock”  was  selected 
as  the  first  alternative. 

BSG  has  been  adopted  as  the  waste  form  to  be  used  at  the  Hanford 
site  [36]  [37].  This  decision  was  based  largely  on  the  analysis  in  the  1982 
report,  which  in  turn  heavily  relied  on  the  recommendations  of  a  peer  review 
panel  (non-DOE  research  scientists)  chaired  by  L.L.  Hench  and  a  report  by 
Schultz  (1980).  It  is  not  certain  that  the  Hench  panel  (or  an  equivalent 
panel  assembled  in  the  late  1980s)  would  have  selected  BSG  for  the  Hanford 
site  waste  (R.  Ewing,  member  of  the  Hench  Panel,  personal  communication, 
1994),  given  that  commercial  scale  production  of  synrock  containing  non¬ 
radioactive  waste  simulants  now  has  been  demonstrated  [38].  An  excellent 
review  of  the  alternative  waste  forms  at  the  time  when  the  Hanford  decision 
was  made  is  [39]. 
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The  acceptance  of  BSG  as  the  sole  waste  immobilization  route  was  jus¬ 
tified  based  on  three  main  assumptions: 

(1)  glasses  accommodate  variations  in  waste  composition 

(2)  glasses  are  reasonably  easy  to  process  and  vitrification  technology  is 
reasonably  well  developed 

(3)  glasses  are  reasonably  durable. 

The  general  validity  of  each  of  the  three  assumptions  noted  above  is 
discussed  briefly  below.  The  BSG  decisions  may  reflect  in  part  the  fact  that 
the  NRC/EPA  regulations  do  not  require  primary  reliance  for  waste  contain¬ 
ment  to  be  placed  on  the  waste  form.  However,  the  long-term  hydrology  and 
stability  of  a  geologic  repository  are  difficult  to  predict,  so  minimization  of  re¬ 
liance  on  repository  design  would  be  preferable.  The  need  for  a  scientifically 
based,  coordinated  program  of  research  to  develop  a  selection  of  second  gen¬ 
eration  alternative  waste  forms  with  enhanced  durability,  minimized  waste 
volume,  lower  exposure  of  humans  to  radiation  during  processing,  and  lower 
processing  costs  will  be  advocated. 

8.2  The  Problem 


At  the  Hanford  site  alone  there  are  approximately  180  tanks  of  waste 
material  awaiting  solidification  and  storage.  Large  volumes  of  material  are 
also  stored  at  the  SRP  (51  tanks  comprising  34  million  gallons  [40])  and 
West  Valley.  Materials  in  these  tanks  have  low  to  high  levels  of  radioactivity 
and  extremely  complex  and  diverse  chemistries.  The  chemical  speciation. 
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oxidation  states  of  the  elements,  and  the  structures  of  the  phases  are  virtually 
unknown.  In  addition,  100+  metric  tons  of  weapons  grade  Pu  (~  50  tons 
U.S.  waste,  50  tons  from  the  former  USSR)  await  immobilization. 

8.3  Borosilicate  Glass  Waste  Loading 


The  amount  of  waste  that  can  be  accommodated  in  a  glass  without 
phase  separation  depends  on  the  nature  of  the  waste  constituents  and  the 
glass  composition  and  structure.  While  glasses  are  “flexible”  networks  (point 
(1)  above)  they  have  quite  Hmited  abilities  to  accommodate  important  sub¬ 
stitutions.  For  example,  the  likely  solubility  of  P,  a  common  constituent  of 
waste  streams,  in  BSG  may  be  less  than  1%.  The  Hanford  waste  contains 
sufficient  phosphorus  that  it  will  have  to  be  removed  chemically.  Otherwise, 
instead  of  the  8,000  logs,  the  tank  waste  will  have  to  be  processed  into  ~ 
200,000  logs  (each  25  cubic  feet)  at  the  cost  of  $1,000,000  per  log  (unoffi¬ 
cial  estimate,  PNL).  Strategies  to  reduce  BSG  volumes  involve  separation 
from  the  sludge  of  constituents  that  degrade  glass  stability  (e.g.,  Cs  and  Sr 
chlorides  and  fluorides  by  ion  exchange:  PNL  is  conducting  a  small  amount 
of  research  in  this  area).  However,  the  method  of  immobilization  of  sepa¬ 
rated  Cs-  and  Sr-rich  materials  is  not  clear  if  BSG  remains  the  only  available 
waste  form.  Unless  effective  separation  strategies  are  developed,  the  cost  and 
attendant  storage  requirements  will  be  of  great  concern. 

The  current  strategy  for  waste  disposal  is  designed  around  manipulation 
of  the  waste  stream  to  suit  the  nature  of  the  BSG  product.  An  alternative 
approach  would  be  to  develop  a  selection  of  waste  forms  suited  to  varia¬ 
tions  in  the  waste  stream.  For  example,  alternative  glass  compositions  (e.g.. 
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Fe-phosphate  glasses,  more  aluminum-rich  glasses)  might  accommodate  un¬ 
modified  or  less  modified  wastes,  reducing  costs  and  human  exposure  during 
processing. 

The  above  discussion  highlights  two  points:  (1)  the  waste  streams  must 
be  carefully  characterized  so  that  strategies  can  be  designed  to  separate  el¬ 
ements  that  have  low  BSG  solubilities  or  undesirable  consequences  for  the 
glass  structure,  and  (2)  alternative  waste  forms  are  needed  to  allow  efficient 
disposal  of  the  full  range  of  waste  materials.  These  research  areas  are  clearly 
important,  and  should  receive  high  priority. 

8.4  Processability 


While  under  ideal  conditions  it  is  likely  that  BSG  could  act  as  a  rela¬ 
tively  effective  barrier  to  radionuclide  release,  the  complexity  of  the  waste 
to  be  stored  and  processing  variables  (e.g.,  BSG  performance  is  greatly  af¬ 
fected  by  flow  rate  [39])  make  this  outcome  improbable.  In  part,  the  initial 
selection  of  BSG  reflected  the  perception  that  it  was  a  “proven”  technology 
(point  (2)).  Processability  was  taken  into  consideration  by  the  Bench  panel 
(1981;  Ewing,  written  communication,  1994)  and  is  discussed  extensively  in 
the  DOE  white  paper  [37].  At  the  time  of  selection  there  was  no  operat¬ 
ing  BSG  plant  in  the  U.S.  The  successful  French  vitrification  operation  at 
La  Hague  was  used  to  support  the  choice  of  BSG.  However,  it  should  be 
noted  that  the  two  operating  lines  at  La  Hague  are  dealing  with  less  com¬ 
plex  commercial  reactor  fuel  under  conditions  where  the  waste  stream  can 
be  readily  manipulated.  Thus,  the  choice  to  use  BSG  for  SRP  and  Hanford 
waste  immobilization  was  not  based  on  experience  with  U.S.  wastes.  Partly 
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for  technical  reasons  associated  with  the  characteristics  of  the  waste  stream, 
the  SRP  vitrification  plant  is  now  bilhons  of  dollars  over  budget  and  5  years 
behind  schedule  [40].  As  of  this  date,  the  SRP  vitrification  facility  has  not 
produced  any  BSG  product.  Given  the  details  noted  in  the  GAO  report, 
concern  over  the  likely  outcome  of  the  Hanford  BSG  operation  is  warranted. 

As  experience  with  storage  of  real  high-level  waste  accumulates,  the 
effectiveness  of  BSG  immobilization  must  be  reexamined.  This  process  of 
reevaluation  appears  to  be  occurring  elsewhere.  For  example,  the  French, 
Chinese,  Japanese,  and  Russians  have  all  expressed  interest  in  synrock  and 
intend  to  use  it  for  wastes  not  suited  to  BSG  (K.  Smith,  ANSTO  Austraha, 
written  communication,  1994).  It  is  of  great  concern  that  the  United  States 
does  not  have  a  viable  research  and  development  program  for  second  gen¬ 
eration  nuclear  waste  forms,  since  a  broader  set  of  options  for  waste  forms 
(including  BSG)  might  well  reduce  costs  and  improve  performance  in  the 
immobilization  of  nuclear  wastes. 


8.5  Durability 


While  durability  is  a  vital  requirement  for  a  nuclear  waste  form  (point 
(3)),  long-term  projections  are  difficult.  Based  on  relatively  short-term  ex¬ 
periments,  numerous  studies  have  argued  that  several  modern  ceramic  waste 
forms  have  much  enhanced  resistance  to  corrosion  compared  to  BSG.  This 
should  not  be  surprising  as  glasses  are,  by  their  very  nature,  metastable  com¬ 
pounds.  Both  comparative  experimental  studies  and  examination  of  natu¬ 
rally  and  experimentally  reacted  natural  analogs  verify  that  glasses  may  be 
easily  altered  at  low  temperatures  over  short  times.  It  should  be  noted  that 
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natural  glasses  do  not  have  boron-rich  compositions  and  that  their  persistence 
in  the  near-surface  environment  reflects  their  limited  exposure  to  water  and 
not  their  inherent  stability.  It  is  widely  known  that  glass  solubility  and  sta¬ 
bility  can  be  degraded  by  various  chemical  substitutions.  Among  the  more 
serious  possible  consequences  of  devitrification  and  phase  separation  is  en¬ 
hanced  solubility  of  waste  constituents  in  solutions  should  groundwater  come 
in  contact  with  the  stored  material.  Reimus  et  al.  [41]  report  40%  increases 
in  leach  rates  due  to  devitrification  of  FY  1987  BSG.  Physical  degradation 
through  cracking  may  also  accompany  crystallization. 

BSG  disposal  strategies  have  been  designed  around  the  inherent  insta- 
bihty  of  BSG  at  high  temperatures  [42]  [43]  [44].  Over  relatively  short  times 
(geologically  speaking),  nuclides  such  as  90  Sr  and  137  Cs  (half  lives  of  28.1 
and  30.2  years  respectively)  generate  heat  within  glasses.  Waste  loadings  of 
glasses  are  severely  limited  by  thermal  instabilities  of  BSG.  This  necessitates 
a  two-stage  storage  process,  places  constraints  on  acceptable  repository  con¬ 
ditions,  and,  inevitably,  contributes  to  higher  costs  due  to  greater  volumes 
of  material  to  be  processed  and  stored. 

8.6  Alternatives  and  Future  Developments 


Based  on  the  work  conducted  up  until  the  time  funding  for  research  on 
alternative  waste  forms  was  discontinued  and  BSG  adopted,  it  is  clear  that 
more  corrosion-resistant,  lower-volume,  cost-effective  alternatives  to  BSG  can 
be  developed  [39].  Although  long-term  nuclear  waste  isolation  strategies 
would  benefit  from  the  availability  of  a  selection  of  waste  forms,  it  seems  likely 
that  BSG  will  remain  the  waste  disposal  method  of  choice  unless  sufficient 
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research  is  done  to  prove  that  alternative  technologies  can  be  developed. 


For  the  future,  an  active  research  effort  to  design  a  selection  of  addi¬ 
tional  waste  forms  based  on  a  thorough  understanding  of  the  nature  of  the 
problem  (composition,  crystal  structure  of  constituent  phases,  particle  size 
(including  colloidal  constituents),  fluid  speciation,  system  phase  relations, 
etc.)  seems  essential.  The  first  indication  that  alternative  waste  forms  are 
being  considered  seriously  comes  from  the  very  recent  initiation  of  a  ceramics 
feasibility  program  at  Lawrence  Livermore  Laboratories  (start  date  June  1, 
1994).  However,  given  the  magnitude  of  the  effort  required  and  the  extraor¬ 
dinarily  large  projected  costs  involved  in  immobilization  of  existing  wastes, 
a  concerted  effort  from  a  team  of  scientists  and  engineers  analogous  to  that 
assembled  for  the  Manhattan  Project  may  be  warranted.  It  is  perhaps  valid 
to  wonder  why  the  new  Environmental  Molecular  Sciences  Laboratory  at  the 
Hanford  site  has  not  been  directed  to  concentrate  significant  effort  in  this 
area. 
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9  ENVIRONMENTAL  AND  MOLECULAR 
SCIENCES  LABORATORY  (EMSL) 


The  EMSL  is  currently  under  construction.  When  completed,  it  will 
provide  Pacific  Northwest  Laboratories  (PNL)  and  the  Department  of  En¬ 
ergy  (DOE)  with  a  magnificently  equipped  laboratory.  The  stated  mission 
of  the  laboratory  is  “to  advance  molecular  science  and  support  the  long¬ 
term  missions  of  the  U.S.  Department  of  Energy.”  As  originally  conceived, 
the  work  in  molecular  science  was  to  be  “far-reaching,  not  directly  linked  to 
environmental  restoration  or  waste  management.” 

The  decision  to  build  EMSL  was  made  at  a  time,  the  mid-1980s,  when 
an  expansive  philosophy  motivated  DOE  managers.  The  Superconducting 
Supercollider  (SSC)  was  proposed  at  this  time,  as  well  as  major  new  facilities 
at  each  of  the  DOE  laboratories.  EMSL  was  to  be  PNL’s  prize. 

As  federal  funding  for  research  becomes  increasingly  tighter,  it  is  im¬ 
perative  to  optimize  every  possible  facet  of  financial  and  human  resources  to 
meet  national  needs.  The  restoration  of  waste  sites  is  a  major  national  pri¬ 
ority  whose  costs  are  unknown,  but  could  be  measured  in  trillions  of  dollars 
for  all  waste  sites.  The  waste  on  DOE  lands  constitutes  a  significant  fraction 
of  the  national  problem.  In  1989,  the  estimated  cost  of  cleaning  up  one  DOE 
site,  Hanford,  was  $40  billion.  Just  three  years  later,  the  estimate  doubled 
to  $80  billion,  with  more  than  30  years  required  to  achieve  the  cleanup. 

The  presentations  with  respect  to  EMSL  reveal  clearly  the  tensions  be¬ 
tween  the  view  that  molecular  science  should  be  pursued  for  its  own  sake,  and 
the  view  that  the  work  of  the  laboratory  should  concentrate  on  the  science 


83 


of  cleanup.  The  magnitude  of  the  cleanup  problem  was  emphasized,  along 
with  the  irrationality  of  some  of  the  regulations.  But  clearly,  the  science 
problems  associated  with  molecular  dynamics,  radiation  chemistry,  surface 
science,  etc.,  are  the  primary  interest  of  the  laboratory. 

A  laboratory  with  an  uncertain  mission  is  bound  to  fail.  Attempting  to 
do  basic  science  as  well  as  cleanup  science  will  inevitably  lead  to  a  situation 
in  which  the  strong  members  of  the  laboratory  are  working  at  the  frontiers 
of  molecular  science,  while  the  lesser  scientists  are  given  the  various  applied 
chores  of  looking  at  such  issues  as  tank  safety  and  subsurface  transport  of 
contaminants. 

We  are  also  concerned  about  the  balance  of  disciplines  within  EMSL.  To 
the  outsider,  the  laboratory  appears  to  have  been  captured  by  the  physical 
chemists.  Physical  chemistry  is  an  important  subject  in  its  own  right,  and 
has  much  to  contribute  to  the  science  of  cleanup.  But  the  biological  sciences 
appear  to  have  been  neglected.  Yet  many  of  the  most  challenging  problems 
of  both  molecular  science  and  cleanup  involve  modern  biology. 

In  these  times  of  budget  stringency,  DOE  cannot  afford  to  support  a 
laboratory  whose  primary  mission  is  to  further  molecular  science  while  sec¬ 
ondarily  working  on  cleanup.  The  mission  of  EMSL  should  be  clearly  stated: 
that  of  developing  the  science  and  technology  of  cleanup.  The  laboratory’s 
performance  and  future  funding  should  be  determined  by  the  progress  the 
laboratory  has  made  toward  developing  workable  technologies  for  cleanup. 
As  indicated  in  our  overview  of  the  Subsurface  Science  Program,  the  scien¬ 
tific  problems  associated  with  cleanup  are  challenging.  They  involve  surface 
science  in  complex  media,  problems  of  genetically  engineering  bacteria,  and 
issues  of  interaction  of  fluids  with  the  subsurface.  These  challenges  should 
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provide  a  situation  in  which  first-class  scientists  work  in  an  atmosphere  with 
good  morale  and  high  scientific  motivation. 


9.1  Recommendations  for  EMSL 


As  described  above,  the  EMSL  is  in  a  difficult  process  of  making  a  tran¬ 
sition  from  non-focused,  open-ended  research  to  mission-oriented  research.  It 
is  clear  that  if  the  transition  is  successful,  the  EMSL  could  make  substantial 
contributions  to  the  remediation  efforts  being  sponsored  by  DOE.  This  is 
because  the  problems  being  encountered  in  the  toxic  waste  cleanup  are  far 
beyond  the  scope  of  solutions  available  from  the  normal  problems  of  indus¬ 
trial  chemical  processing.  If  the  cleanup  is  to  be  successful  at  a  cost  that 
does  not  bankrupt  the  nation,  innovative  solutions  must  be  developed.  The 
role  of  basic  research  in  the  process  should  be  to  consider  a  broad  range  of 
novel  possible  processes  relevant  to  cleanup.  For  this  process  to  be  useful,  the 
ideas  developed  in  basic  research  must  be  transferred  to  a  practical  testing 
stage  and  evaluated  for  real  utility.  If  truly  novel  and  innovative  ideas  are 
pursued,  it  is  clear  that  only  a  small  fraction  of  the  ideas  will  be  practical, 
and  evaluation  of  the  success  of  the  basic  research  effort  must  take  this  into 
account.  However,  it  is  also  clear  that  successful  application  of  the  ideas 
developed  in  basic  research  can  only  occur  if  the  basic  research  scientists  are 
responsible  for  communicating  their  ideas  to  those  involved  in  evaluation, 
and  seeing  the  ideas  through  the  evaluation  process. 

The  organizational  charts  for  the  EMSL  and  its  function  within  PNL  are 
shown  in  Figures  8-1  and  8-2.  The  activities  of  EMSL  are  organized  as  four 
fundamental  research  sub-groups.  Theory  Modeling  and  Simulation,  Chemi- 
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Figure  9-1 .  A  multidisciplinaty  laboratory  for  Environmental  Molecular  Science. 
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Advanced  Processing  Characterization 


cal  Structures  and  Dynamics,  Macromolecular  Structure  and  Dynamics,  and 
Materials  and  Interfaces,  which  are  to  transfer  basic  research  results  to  two 
bridging  programs.  Environmental  Dynamics  and  Simulation  and  Advanced 
Processing.  The  two  bridging  programs  are  expected  to  develop  apphca- 
tions  of  the  basic  research  results  that  can  be  used  in  the  remediation  efforts 
being  carried  out  in  other  laboratories  in  PNL.  As  noted  above,  the  organi¬ 
zational  plan  places  an  inadequate  emphasis  on  microbiology,  which  has  the 
potential  for  significantly  improving  remediation  capabilities.  To  encourage 
necessary  interdisciplinary  work,  microbiology  should  be  added  within  the 
existing  groups  (perhaps  with  modifications  of  their  present  titles),  rather 
than  as  a  separate  group.  In  addition,  the  EMSL  does  not  yet  appear  to  be 
functioning  fully  as  envisioned  in  the  present  organizational  charts,  and  we 
present  below  some  suggestions  as  to  how  operations  could  be  improved.  Suc¬ 
cessful  implementation  of  these  suggestions  will  require  managerial  initiative 
to  foster  teamwork  and  coordinate  research  directions. 


9.2  Advanced  Processing 


The  Advanced  Processing  Group  is  primarily  tasked  with  developing 
apphcations  for  the  cleanup  effort  of  tank  wastes.  In  this  problem,  all  the 
normal  issues  of  chemical  processing,  such  as  materials  transport,  separation, 
and  selective  reaction  are  complicated  by  the  nature  of  the  mixed  wastes 
that  must  be  processed.  The  Advanced  Processing  Group  appears  to  be  well 
placed  in  awareness  of  the  real  problems  involved  in  tank  waste  processing, 
and  of  the  nature  of  the  approaches  that  might  be  useful  in  dealing  with  the 
mixed  wastes.  However,  the  amount  of  input  that  it  is  receiving  from  the 
basic  research  groups  on  these  problems  appears  to  be  limited.  In  some  of  the 
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sections  below,  we  will  comment  on  how  this  could  be  improved  by  efforts  of 
the  basic  research  groups.  One  serious  omission  in  the  efforts  of  the  Advanced 
Processing  Group  is  that  there  appears  to  be  little  or  no  effort  in  evaluating 
microbiological  processes  for  toxic  waste  degradation.  Such  processes  have 
been  developed  successfully  in  industrial  apphcations,  and  should  be  tested 
for  the  difficult  problems  that  must  be  addressed  by  PNL. 

9.3  Environmental  Dynamics  and  Simulation 

The  Environmental  Dynamics  and  Simulation  Group  has  developed  a 
well  focused  program  using  the  advanced  computing  capabilities  of  the  model 
reactive  chemical  transport  in  contaminated  soils  and  ground  water.  The 
group  has  developed  active  collaborations  with  the  Subsurface  Program,  and 
should  in  time  broaden  the  scope  of  its  work  to  include  bacterial  transport. 
There  does  not  yet  seem  to  be  much  input  from  the  basic  research  efforts 
within  EMSL,  which,  as  noted  in  earlier  sections  of  this  report  and  below, 
will  be  needed  in  developing  simulations  that  provide  realistic  assessments  of 
what  is  happening  in  the  variable  subsurface  environment. 

9.4  Theory,  Modeling  and  Simulation 

The  Theory,  Modeling  and  Simulation  Group  has  developed  capabilities 
to  realize  the  long-held  dream  of  molecular  modeling  to  tailor  materials  to 
desired  functions.  The  accomplishments  of  this  group  in  developing  the  abil¬ 
ity  to  calculate  structure  and  function  for  complex  molecules  such  as  those  in 
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zeolites  are  impressive.  However,  without  follow-through  to  see  that  the  cal¬ 
culations  are  tested  experimentally,  the  work  of  this  group  cannot  be  useful 
to  the  overall  mission  of  PNL.  It  is  extremely  surprising  that  there  appears 
to  be  little  or  no  correlated  experimental  work  within  the  three  experimental 
basic  research  groups  to  test  the  predictions  of  the  calculations.  The  success 
of  EMSL  requires  that  systematic  collaborations  of  theory  and  experiment  be 
established  to  predict  and  test  the  properties  of  materials  that  are  relevant 
to  toxic  waste  separations  and  degradation.  In  this  regard,  the  calculational 
efforts  in  rational  enzyme  design  could  have  a  significant  positive  impact  in 
helping  to  develop  microbiological  processes  for  tank  waste  remediation,  if  it 
were  coupled  to  correlated  experimental  studies.  This  group  needs  to  place 
“Developing  Real  Collaborations  with  Experimentalists”  at  the  top  of  its  Ust 
of  “Challenges.” 


9.5  Materials  and  Interfaces 


Materials  and  interfaces  axe  clearly  of  great  importance  in  both  the 
separations  and  catalysis  required  for  tank  remediation,  and  in  the  transport 
processes  that  must  be  understood  for  cleaning  up  ground  contamination. 
The  Materials  and  Interfaces  Group  has  developed  experimental  capabilities 
that  could  be  applied  in  both  of  these  areas,  and  has  begun  to  develop  an 
effort  focused  on  the  mission  of  PNL.  Immediate  opportunities  are  open  for 
such  a  focused  eflFort  in  at  least  two  areas.  One  would  be  in  synthesis  and 
characterization  of  the  materials,  such  as  zeolites  and  other  catalysts,  being 
studied  by  the  Theory,  Modeling  and  Simulation  Group.  Another  would  be 
in  investigating  the  nature  of  adsorption  and  adhesion  on  the  highly  variable 
surface  structures  of  the  subsurface  environment.  The  latter  area  of  research 
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should  include  microbiological  investigations  as  well  as  chemical  studies.  This 
work  would  be  usefully  addressed  with  scanned  probe  microscopies,  as  these 
techniques  require  neither  a  vacuum  environment  nor  spatial  homogeneity  of 
the  sample.  Both  of  these  areas  of  research  require  serious  collaboration  with 
the  prospective  users  of  the  information  generated  in  order  to  be  useful. 

9.6  Chemical  Structures  and  Dynamics 


The  Chemical  Structure  and  Dynamics  Group  has  developed  state-of- 
the  art  capabilities  in  investigating  reaction  mechanisms  and  kinetics.  The 
initial  demonstrations  of  these  capabilities  have  been  performed  on  subjects 
that  have  general  significance  for  remediation,  but  little  or  no  direct  con¬ 
nection  has  been  made  to  other  activities  in  EMSL  or  PNL  even  where  the 
applications  seem  obvious,  as  with  radiolytic  reaction  mechanisms.  More 
rapid  progress  toward  making  such  connections  could  be  made  by  collabo¬ 
rating  with  the  Theory  and  Modeling  Group  to  test  predictions  made  for 
both  solid  state  and  enzyme  catalysts,  and  then  to  work  with  the  Advanced 
Processing  Group  in  transferring  the  most  promising  materials  to  the  testing 
stage.  Another  area  in  which  the  capabilities  of  this  laboratory  could  be 
applied  is  in  state-of-the-art  characterization  of  the  structure  and  function 
of  bacteria  that  may  be  used  in  bioremediation.  It  would  be  natural  and 
effective  to  add  microbiologists  to  the  staff  of  the  Chemical  Structure  and 
Dynamics  Group  to  encourage  such  interdisciplinary  work. 
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9.7  Macromolecular  Structure  and  Dynamics 


The  Macromolecular  Structure  and  Dynamics  Group  has  developed  and 
is  continuing  to  develop  state-of-the-art  NMR  capabihties  for  studying  cel¬ 
lular  and  structural  biology.  They  have  developed  collaborative  interactions 
outside  of  EMSL  involving  assessment  of  health  risks  of  exposure  to  toxic 
chemical  and  radiation.  Similar  collaborative  efforts  within  EMSL  should  be 
directed  at  studies  of  enzyme-based  degradation  of  toxic  waste,  and  studies 
of  structure  and  function  of  bacteria  that  might  be  used  in  bioremediation. 
The  studies  of  enzyme  function  should  be  coordinated  with  the  theoretical 
work  of  the  Theory  and  Modeling  Group,  and  with  experimental  studies  in 
the  other  research  groups  of  EMSL.  They  should  be  directed  at  transferring 
results  to  the  Advanced  Processing  Group.  For  the  studies  of  microbiology, 
it  would  be  natural  and  effective  to  add  microbiologists  to  the  staff  of  the 
Macromolecular  Structure  and  Dynamics  Group,  perhaps  with  joint  mem¬ 
bership  in  the  Chemical  Structure  and  Dynamics  Group. 
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10  BIOREMEDIATION 

10.1  The  Long  Haul 


The  program  falls  naturally,  it  appears,  into  at  least  two  parts — those 
dealing  with  short-term  scientific  problems  and  those  dealing  with  longer- 
term  ones.  The  horizon  for  the  short-term  problems  is  optimistically  set  at 
about  ten  years.  The  longer  ones  may  take  virtually  forever. 

There  is  a  substantial  probability  that  there  are  genuine  short-term  sci¬ 
entific  problems — that  the  distribution  and  heterogeneity  of  soil  and  subsoil 
types  is  terribly  intricate  and  poorly  structured  taxonomically,  that  soil  hy¬ 
drology  is  terribly  complex,  that  the  transport  of  bacteria  through  subsoils  is 
too  varied  a  problem  to  effectively  deal  with,  etc.  Then  there  are  also  all  the 
problems  associated  with  the  bacteria  themselves:  their  metabolic  activity 
in  extreme  environments,  their  tendency  to  adhere,  etc. 

These  are  the  “short-term”  problems.  The  long-term  problems  address 
the  origins  of  life  itself. 

If  there  are  no  genuine  short-term  problems,  where  will  the  program 
find  itself  in  10  to  20  years?  A  great  deal  of  scientific  research  will  have 
been  done,  much  of  it,  it  is  true,  of  broad  interest  and  importance.  But 
no  scientific  problem  of  central  relevance  to  the  problem  of  bioremediation 
will  have  been  dented,  much  less  solved,  and  certainly  the  whole  group  of 
problems  that  must  be  solved  collectively  according  to  the  organizational 
goals  of  the  program  may  hardly  have  been  touched. 
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Is  the  program  to  be  viewed,  then,  as  a  very  long-term  research  under¬ 
taking  which  at  some  unspecified  time  far  in  the  future  may  solve  some  of 
our  pollution  problems?  Or  are  there  relatively  short-term  and  reasonable 
goals  for  the  program  by  means  of  which  some  positive  progress  can  be  seen 
along  the  way? 

Field  tests  at  near  or  full  scale  seem  appropriate.  Useful  experimental 
data  will  be  of  great  significance  to  the  long-term  program  not  only  as  a 
guide  along  the  way  to  validate  predictions,  but  also  as  essential  ingredients 
of  the  actual  efforts  at  bioremediation. 

Normally,  near  full-scale  tests  are  inappropriate  to  a  subject  just  barely 
in  its  infancy.  But  in  this  particular  instance,  the  program  must  early  on  face 
the  questions  of  licensing  and  making  statements  of  environmental  impact. 
If  these  cannot  be  overcome  now  at  a  purely  harmless  experimental  level, 
how  can  they  possibly  be  met  when  it  is  proposed  to  introduce  living  organ¬ 
isms  in  bulk  into  deep  terrain,  possibly  organisms  that  have  been  genetically 
modified,  and  whose  full  range  of  possible  appetites  and  mutants  is  hardly 
known — indeed  it  is  never  known.  Moreover,  the  organisms  will  be  operating 
far  out  of  view  and  oversight. 

The  obvious  large-scale  experiments  address,  for  example,  questions  of 
hydrology,  diffusion  and  behavior  of  acquifers.  We  would  suggest  a  large-scale 
and  patient  examination  of  diffusion  of  water  in  typical  terrains,  especially 
the  types  of  terrains  where  remediation  is  obviously  called  for  and  perhaps 
even  at  the  sites  themselves.  This  will  involve  pumping  water  with  tracers 
into  the  ground,  and  watching  over  a  long  period  of  time  its  arrival  at  loca¬ 
tions  at  several  scales  of  remoteness.  This  kind  of  experiment  could  probably 
be  carried  out  at  Hanford  without  upsetting  anyone,  and  the  data  collected 
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will  be  of  immense  value.  No  experiment  with  glass  beads  in  a  laboratory  set¬ 
ting  will  reveal  the  deep  hydrology  of  a  particular  site  as  well  as  experiments 
at  that  site  itself. 

If  it  should  come  to  pass  that  hydrology  at  some  polluted  site  is  com¬ 
fortably  understood,  then  appropriate  nutrients  and  relatively  benign  {but 
effective)  ordinary  organisms  might  perhaps  be  introduced  without  much  le¬ 
gal  difficulty,  though  we  wonder  about  this.  One  cannot  help  but  recall  the 
raging  controversy  about  disposal  of  radioactive  wastes  in  salt  domes,  though 
the  temporal  scale  of  this  proposal  is  of  altogether  different  magnitude,  it  is 
true. 


Let  us  summarize,  for  emphasis,  what  we  are  driving  at  above.  The 
bioremediation  program  is  one  of  tasked  pure  research,  or  pure  research  with 
a  mission.  The  general  wisdom  of  science  is,  however,  that  such  an  effort 
cannot  be  managed  on  a  schedule.  Significant  scientific  results  emerge  in 
their  own  good  time. 

To  give  the  bioremediation  program  a  better  transitional  form,  we  think 
it  desirable  that  large-scale  field  experiments,  in  anticipation  of  remedial 
progress  and  insights,  be  conducted  with  all  vigor  along  the  way. 
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